General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




UNIVERSITY OF 




SOUTHERN CALIFOEINIA 


CIVIL ENGINEERING DEPARTMENT 


FINITE ELEMENT ANALYSIS 


OF LAMINATED PLATES AND SHELLS 

(NASA-CE-157107) FINITE ELEMENT ANALYSIS CF N78-24578 

LAMINATED PLSTES AND SHELLS, VOLUME 2 
Annual Progress Eeport, Mar. 1977 - Mar. 

1978 (University of Southern Calif.) 164 p Unclas 

HC A08/MF A01 CSCL 20K G3/39 20793 


PETER N.-H. CHANG 


ANNUAL PROGRESS REPORT 


MARCH 1977 - MARCH 1978 




Volume 2 


NASA Grant No, 1401 



University of Southern California 
University Park 
Los Angeles, California 90007 



FINITE ELEMENT ANALYSIS OF 
LAMINATED PLATES AND SHELLS 


Paul Seide 
and 

Peter N.-H. Chang 


Annual Progress Report 
March 1977 - March 1978 


Volume 2 



NASA Grant No, NSG 1401 


t lw i j ) F.n » - 3 L im ^ n ^ r ^ rtiTiiii^"^ tjihii E i i ill ^ t ? - -' &BT- . zrT^ igu.:. i i i ■‘■ ■ ^ 


APPEiroiX A* COI-ffUTER HOCEAM FCE lAIOTM AMD CIUOTICAL SHEII, AB-LKSIS 


L.ASflE^ aFiltri 


FiJRTRAM V*5At6Sl) /Ki 2:>-?flAif-7d 


1 2£^u PAut^ 1. 


00010 
Q0020 
00030 
00040 
OOC50 
00060 
00070 
OOOBO 
00090 
OOlOO 
001X0 
00 120 
00130 
00140 
00150 
00160 
00170 
OC 160 
00190 
00200 
00210 
00220 
00230 
00240 
00250 
00260 
00270 
00200 
00300 
0Q310 
00320 
Q0330 
0034C 
00350 
00360 
00370 
00360 
00390 
00400 
00410 
00420 
00430 
QU440 
0045C 
00460 
0047 0 
00400 
00490 
00500 
00510 
00520 
00530 
Q0 540 
00550 
00560 

ousro 


lMPulC.1T Ru4L*0( A-H.u-Z) 

REAU TtTT.S.3T.HE0 

LOaiCAu £UPRT*EuPCH*ufcNPRT«GEli^CH,w»sFP^rtt6uuST 

LOMMQN V6AR/ aC 99991 

CUHHQW /30NA /libOllb] ,JUK1400» 

COHHaN /EuPA'R/ MPAR ( 1 4 1 • NUHNP *46 ANJ. NE..T TP *N 1 * M2* H3 * N4 >n5 * HT6T , MEO 
CQHHQN /EM / 000(25471 

COHMCM /OYN/ roUallil *NUVU*Trt6Purt»Tri63Crl 
COMMON /TAPES/ MQOlol 

CCHMON /EXTRA/ HOOEA.N t d*N lO « NT 1 ^ iuEu3 «NUMEU » r ( 12 ) «5 (12 1 
COMMON /30C/ NQLOCK,*NEQit>*uuiMF* l./3HtNslu*NAa«NVVtAN0RM*Nfja 
common /CTL/ 0RAV»eLPHT.EuPCH*3£NPKT*MuNP2N.UEFr*Crt»lJilTC6. 

J MINONP* IPET iNMuUTa tNP*-T . JNAAl lOl 

COMMON /o 6U5TF/ 6EiJST«MB.E6EO 
C0HMaN/SElS/NRSCtKAP6.NRS«IT*0JATE*JFUA6 


OPENlUNlT- 
□PEN(UN1T' 
OPEN(UNlT- 
□PENIUNIT: 
□PEN (UNIT- 
OPEN (UNIT- 
QPEN(UN1T= 
OPEN (UNI TS 
OPEN (UNITS 
□PErt(UNlT= 
OPEN (UNI TS 
OPEN (UNI Is 
PPEN< UNITS 
UPENIUNITs 


1 ,DEV ICE= 
2 » DEVICES 
3f0dV ICEs 
4. DEVICES 
5»UEVECEs 
6. DEVICES 

7 .DEVICES 

6 a DEV 1 CEs 

9 .DEVICES 
to. DEV ICES 
12 » DEV ICEs 

lo * DEVICES 

14*D£VlCus 
15 .DEVICES 


»DSK» 

»USK* 

*D3M* 

*D3K* 

►DbA* 

»DSK" 

•05K-* 

»D5K* 

•DSK* 

*USH* 

• DSK» 
» J3K* 
•USA* 

• D3K* 


•FluEs 
*Fli>Es 
*F I UCs 
♦FILLS 
.FiuEs 
t fills 
♦ FIuEs 
iF I uL= 
»F ILus 
.FILES 

.FILLS 

.files 

*F 1 LLs 
.FILES 


•CCOl 

•CL02 

• CC03 
■cca4 
«a205 

• 3206 
*uC07 
•CL33 
•CC09 
•CClO 
«CC12 
•CC14 
*Ctll4 
•CCI5 


.DAT * 1 
.DAT •) 
.DAT • ) 
*DrtT* 1 
. DAT * 1 

• DAT*) 
. DAT 0 1 

• DAT*) 

• DAI * ) 
•DAT*} 

• DAT*) 

• DaT ’ ) 
.DAT " ) 
.DAI ^ > 


NPLOT3=0 

a CALL HEUPRT(KUVN) 
JFLAiJ=0 
■RITE (13) HED 


INPUT 


U U 1 N r 


DATA 


N2s:N1'<-64NUHNP 

N3=N2FNUMNP 

N4SSN3-I-NUMNP 

N5=W4+NUMNP 

N0=:N5FNUMNP 

CALL ERKUR(N61 

read Cb.lCOll NPAR 

uALu InPJTJ (A(N1 ) iA(N2) rA(Na). a(N 4I *A( Mb) 1 
STORE THE NODAL CODRIUINATLs FJR L.\TER JbE 


2 CZT] ^3 


3 ^ 


e^ , E3 


Vj 


I 


MAIN. 

LASHEL.FtJM FOrtTKAN V.5A(621) /f.1 33“«A1f-7B 

t2;oa PAGE 

> oosao 


WRITE 11;}) t A[ n .l^NEf N5) 


00590 

C 



00600 

C 

FDRM E1.CMENT 3T1FFMC.5 8LS 


0061C 

c 



) 00620 

c 



00630 


MtiANi3=0 


C0640 


NUMBt_=0 


0065 0 

■ c 



1 00660 


jATA PQRTHCJLb SAVE 


00670 


IFIHdOEX.EQ.il WHITE WTa»4J NPAH 


00680 


NT 1=1 


00690 


If IMlNBNO .8E- 1 .ANOi MOOEX .EQ. 0) NT 1=4 


) 00700 


IF {HJNEnO *EE* 1 »AN0* HUOEa « {^3 . 1.J i.T 1= 7 


00710 


WRITE (NTlI NPAR 


00720 


NUHEE=NdHEL4NPAR {21 


00730 


MTYPE=NPAft£l I 


1 00740 

c 



00750 


CALi- ELTYPEtrtTVPEJ 


00760 

c 



00770 

c 



’ 00780 

c 

ADO an ENO-QF-FILE F(-A0 ON F Il.5 13 


00790 

c 



00000 


H=-l 


ooaiu 


IF tiPET -ET. 0) 


• 00820 


XWRXTEC131 H.M 


00030 

c 



00040 

c 

. 


00650 

c 

DETERMINE duOCAalZE 


00060 

■ c 



00070 

- c 

ADU5TF 


oooao 

c 



00090 


NEOO=(HTUT “ 4#EE)/lHdAN0 + i-i_ + 11/2 


00900 

c 



cosio 

c 

DVER-KlUE THE SYSTEM MATRIX ouOEA51ii:.£ wITH THE INPUT iN0N*-2hR0. 

C092O 

c 

VAEUE. KEQQ. 


00930 

c 

THIS over-ride ENTRY IS TO ALEJw PRUQRA'3 CHEEKING 

OF MULT I- 

1 00940 

c 

tjECCK aegjkithhs with whai *ooi-o normae».y he one 1 

HLQCK DATA. 

00950 

c 



00960 


if t KEQU.eT.NEgOI NECiH = XdOu 


00970 

c 



009U0 

c 



C099O 

c 

STATIC SDeUTION 


0 1000 

c 



01010 

o90 

CONTINUE 


• 01030 


NEOBl=lHruT - MdANOJ/12’^(MdANOl-E._l + l» 


01030 


NeQu2=IMTUT - H hand - EE» 1 H3ANJ~2 I 1/ 1 3*eL + NOANO 

+ 1 * 

01040 


IF iNEQoi.ET.NEQBl N£QH=NEUl»1 


01050 


IF INEOda.ET.NHOEi NEQB=NE0a2 


1 01060 


NHLUGK = (NeQ-lj/N£ao tl 


01070 


IF (NEOH.uI . nEQI N£Qb=NEU 


01080 


IFtlPET.EQ.OlGLJ TO 693 


01090 


NT dl=204NLJrtNP+NiiQH*NaEUi-K*EL 


j 01100 


IF tMTOl.^jT .NTOriOU TO 790 


0111 0 


WHI TE16. 2020 TNT UT 






IHZE; 


r 


^3 






r " n nzj [zn czD [zz] czz cm cm . cm 


MAIN. 


LASHfcu*Pt]R 


f=UHTHAN V. DA (6^1) VM ^D-MnY-Td 


1^; Jb 


PAbc 1—2 


Q1 lAO 
Cl 150 
01160 
oiiro 

01 ISO 
OllSrO 
012Q0 
□ I Zl 0 
01220 
01230 
012*0 
01250 
01 260 
01270 
01260 
01290 
01300 
01310 
01320 
01330 
01340 
01350 
01360 
01370 
01380 
01390 
0140 0 
01410 
01420 


792 

C 

C, 

c 


lF (MT0T*oT .NTDTJGCj TO 790 
!«RlTt(6,20i0 JNTCT 
MDDEX=l 
GU TO 790 


790 IF (N£Q6 .GE* 2J GO TO 7V2 
IF tNEQ .t-T. 3J GU TO 792 
MRITE (6«2010) NEQO 
NtUtl=2 

^8LULKi:(NcQ-l )/N£mO + I 
HU0EA=1 
CONTTNOE 


1 N PUT 


N U D A L. 


L O A O S 


nse°neqd*nbluck; 

N2=Nl+6*NUHNP 

N3i=N2+LL4NSJ 

NA^Nj^NSd 

call, error 1n*» 

WRITE l6f2Ul) NcQ*MpAiNtl,NE0lj*r<3|_ClCA 

CALL INLlAt Nl ) »AtN21 i A( N3 > * NuMMP tNcOu* SSd «LL j N dL OCAl 

CHElG fur unrestraineu OEGREES QF FREEuOH 

TNEO*NGi'i6LyNci.TfP) 

S T I F r i>t £ o » 


01430 



ND=H2«^REa 

01440 



CALL CHECK! A (Nil 

01450 

C 



01460 

C 


FLK;( TtiTA 

0X470 

c 



01480 



N62u=2*NEUb 

01490 



N2=Nl+NEUd*MdANU 

01500 



N3=tN2+NE0tJit‘L.L 

01510 



N4-N3+44LU 

01520 



NN2=N1+NE26(-H6AN0 

01530 



NN3= NN 2+ NE 2 84Li_ 

01540 



NN4=NN3+4+i.u 

01550 

c 



OlSbO 



CALL AUUSTF (aIN 

01 570 


XU jANURHiNVVJ 

OISOO 

c 



01590 


18 

LUNIlNUh 

01600 

c 



01610 

c 


6 L L U T 1 U N 

01620 

c 



01630 

L 



01640 

L 


STATIC Solution 

01650 

C 



01660 


30 

IFlHUJhX.eUiO ) uU 

0167C 



GO TO 90 

0I6B0 

c 



01690 


32 

RLW INO 12 


PHASE 


s § 

5=! 

i| 

O 
c?d 


cc 


ci 

P. 


l i 




< oiSf’i.r,‘ ** !i - iV/ 'v- •' 


■ I'lBi M~l 'h I l’ II 1l« r I ■■■li/'il - itirV.' 









MAIN. LASttki_.eijH 


fdHTKAN V.sAluaU /*tl a^-rtAVVel 


lajati HaaU 1-J 



01700 

01710 

oi7ao 

01730 

01740 

01750 

01760 

01770 

017B0 

01790 

oiaoo 

01310 

D1320 

oiaao 
01540 
01550 
01B60 
01870 
01880 
01890 
01900 
01910 
01920 
01930 
01940 
019S0 
01960 
01970 
01980 
01990 
02000 
020 1 0 
02020 
O3O30 
02040 
02050 
02060 
02070 
02080 
C209D 
02100 
02110 


CALL 6Ui-fc.Q 

IF i.NOT. GEDbrl iiU TO 90 
WRITE (6*20001 NJMEL*NEeGEIJ 
C . 

C ■ 

<1 . ■ 

c 

90 GO TO 5 

c ■ 

20 I FURWAT (A8mE QUATION PARAMETEKS*// 

X TOTAL NUMBER OF EQUAnaNb =,15* 

X /34H aANOWliiTri =,I6, 

X /34rt NUMbER OF EilUATlUNS IN A bu^LX =.15, 

X /34rt NUMBER OF Di-OCKS =,i5> 

203 FORMAT 11H1,31H0 VERALL TIME LUG*// 

X bXibOHNJOAu POINT INPUT =» F8,2 / 

X 5X*30HELEMc.NT STIFFNESS FORMATION =, F3*2 / 

X SX»30HNOOAL load input =• F8.2 / 

X 5X,30HTOTAL STIFFNESS FORMATION = , Ffl ,2 J' 

X 5X,30rtSTATIC ANALYSIS =* F8*2 / 

X 5X, 30 riel gen VALUE EXTRACTION =t F8.a / 

X SX,30iiSTRUCTURE PLOTTING =. F8 ,2 / 

X bX.jOHFuRCED RESPONSE ANALYSIS =# Fb.ji / 

X SX»3UH«£SP0NSE spectrum analysis =* F8,c ✓ 

X 5X,3QHSTEP- BY-STEP INTEGRATION =, PS .2- / 

X SXiSOrtFREOUENCY RESPONSE ANALYSIS =, F3,2 // 

X SX.oOHTuTAL solution TIME ' =. F8.2 /t 

C 

c ■ . . 

c 

1001 FORMAT (14 la I 

2000 FORMAT ( * OTOTaU NUMaER OF ELEMENTS 

X * NUMBER OF ELEMENTS WITH UEOM^rRi*. ST 1FFNE5S=« , ItiJ 

2010 FUHHAT warning the available bL.ANX CUMMDNT (MTOTl la* * 

A* TOCi SHALL* */,« COMPUTED NO, OF EQUATIONS PER BLOCK- *, 1 10 */, 

X* PROwEEJ IN THE DATA CHECK HDOE WITH NEau=£* i 
2020 FORMAT t*0*44 FATAL ERROR +** AVAlLAbLE CORE Ivs TOO SMALL* 

I,* FOR PuGTTlNG, REnUlREU Ml 07= *.110/ 

2* continue" IN THE OATA CHECK MODE* I 
2830 FORMATCO*** FATAL ERHUR W4W AVAILABLE CORE IS IHAJAQUATE* 
l,» nCOOlKEO MTUT=* • 110/*C0NTINUL IN TriE OATa CHECK MODE*} 

END 


' COMMON BLOCKS 

/GAR/t+47836) 
A +0 

I /JUNK/ (+642) 



HEO +0 

JUK 

+22 






J 

/ELPAR/'i+30) 
NPAR +0 

N2 +22 

NEC *zr 

/EH/ 1+13076) 

NO MNP 
N3 

+ 16 
+23 

MdANO 

N4 

+ 17 
+24 

NElTVP 

NS 

+,.0 

+2S 

N1 

MTPT 


) 




r- 



\kzj'' c~j lzzj cm cm cm ^ ^ ^ cm Em cm Em es 

^ .. 

f ■ ■ . . ■ ■ . 

p , ' ■ , . . ■ ■ 

E ■ ' • ■ ' 



HAiN. LASHEu* 

FUR 

FUKTHAfM 

\y ,541621 1 /Al 

20-.MAY 

-7 0 

1213B 

PAGE 

‘ 

OUQ *0 
/OYN/ 1+201 
I Dub +0 

NuY« 

+13 

Trio PC H 

+14 

I rlSPCrt 

+10 



t 

/TAPES/ (+6) 
NQQ +0 

/EXTHA/1+361 
HDDEX +0 

NTS 

+1 

NIOSV 

+2 

NTia 

+3 

KhOB 

+4 


NUHEL +5 

T 

+6 

S 

+22 






/SOL/1 + ldi> 
NBLQCK +0 

NEvJU 

+1 

LL 

+2 

Nr 

+3 

lUUM 

t4 

1 

NEIG +5 

NAD 

+0 

NVV 

+ 7 

anor>i 

+1C 

NFD 

+ 12 


/CTL/l+2b) 
GRAV +0 

ELPBT 

+2 

ELPLri 

+3 

OENPRI 

*4 

bENPCH 

+5 

» 

OEFPCH +6 

IWTCG 

+7 

HINBNU 

+ 10 

IPLT 

+11 

NPLQTS 

+12 


NPLT +13 

JNKK 

+14 








/GEQSTF/l+Zl 









1 

cEosr +0 

NELGEU 

+i 








/SE IS/ 1+71 
NRSC +0 

KAPG 

+1 

MRS 

+2 

IT 

+3 

DDATE 

+4 

) 

JFLAG . +b 










SUDPKOGKAKS CALLED 








i 

ADDSTF SULEQ 
INPUTJ INL 

HEDPRT 

ELTYPc 

check 

ERRLirl 






SCALARS AND AkKaVS i NO EXPLICIT DEFINJTIuN 

— H 41 I 

iuDT rtc.F 

^RENLED • 1 


J 

• NN2 1 

♦ NTOT 

2 

ASF 


*n 6 

3 

4NN3 

4 


+NE20 5 


l> 

♦ NN4 

7 

«nsd 

10 

•SOOUO 

11 


4NEat32 12 

*NTt 

13 

♦kdyn 

14 

OrtTVPi 

13 

4NEQD1 

lo 

1 

♦ I 17 

XTT 









TEMPORARIES 









i 

■UOOOO 421 

•00001 

422 

•00002 

423 

•U00J3 424 

■U0004 

425 


MAIN, i NL» tRFtUKS JCTc^LfEP J. 
I 




; 


t 

} ■' 




i.A.5HtLol-Ort 


^UHIhtAN VobAtuiill /K1 


;dJ-MAV -7 ct 


oai20 
02S30 
021A0 
02150 
02160 
02170 
02180 
02190 
02200 
02210 
02220 
022 JO 
022AO 
02250 
02260 
0227C 
02280 
0229 0 
02300 
02310 
02320 
02330 
02340 
02350 
02360 
02370 
0 2380 
02390 
02400 
02410 
02420 
C2430 
02440 
02480 
024b0 
02470 
02480 
02490 
0280Q 
02S10 
02520 
0253 0 
02840 
02SS0 
02560 
02570 
02580 
C2S90 
02600 
02610 
C2620 
02630 
02640 
02650 
02660 
C267D 


5UbK6UTlNt: rtDDaTF trt*13fbTrtiTMA55iJ^JHfcl.f>idl-0Clfi»Ni:;2B*l_L»HbANJtANiJHM» 
INVVJ 

IHPI-ICIT M 5 AL 4 dCA--H»U- 2 ) 

' LOGICAL ci-PRT,eLPCri.oeNPRTiGENPCri»OEFPi-H 
RtAL PCT 

CALLEO or MAIN 

FORMS GLLIdAu EQUlLloRIUM EQUATIONS IN duOCKS 

pi HbNSlJU A(Nh 2 n»MaANP} ( 0 tNb 2 uiLuJsSTHl 4 ILL) «TMA 5 StN 8 . 2 bl 

COMMON /OYN/ NTrNur ,A lFA>PTbJETA.NFH.NGM*NAT.NDYN 
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NUHE-NUMEl 
NUM7 =0 
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13U lUll >L}=UllltLl+SSII-fAiv)4iSTH(«j,L) 

DO 150 J=ifNU 
3J=Lm Jl+LHN 
IFIJJ) 1S0.160«132 
132 IF (J-i) 140»|34 i 134 
l.»4 KK=l 1-1 14ND-X41 l-U/2<-J 
Gu ro 145 

140 KK^l J-1 J >I«NU-J4( J-1 1/2+ I 
145 A( I l*JU)=At 1 1*3J 1 + SSIKM 

166 continue 

loO CUNl INUE 



ADUSTF LAiMcL-.J-QR FjKTKAM V.‘£>A(o<in /Kl ld.i’^a HAuu 


032A0 C 

03250 C JEtaHHi.Nii IF 5>TifFNt.Si» IS TQ tJc PLACED JN TAPt: T 

OSZbC C 

05H7Q IF tHM.ifT.i; GO TO 700 

t 03200 DO 200 I-li«4D 

03290 II=LWtn -WSHIFT 

03300 XF{II*GT«NE2S>AN0«11«£-E>ME3DJ GU to 210 

03310 20C CONTINUE 

' 03320 GO 10 700 

03330 21U «niT£(7] LROf NO • [LM( 1) » 1=1 >ND) . (SS t i I a 1^ 1« i.Kcf J 

C334C N0H7=NUM7+1 

Q33S0 t 

03360 700 LJNTlNJii 

03370 UD 710 1.s1«NEC13 

033B0 ANuRM=ANGRM + AiLfU 

03390 IF (A(l» 11 -NE« C.ODO *AND« AlLalJ aLT* ANIN} AH1N=A[|-«U 

03400 IF lAlutU «GT. AHAaI AMAX=A1L»1) 

03410 IF lAlLi U •NC*0 >D0) Nu£G=MDEGF 1 

03420 IF lAlLf 11.EQ.0*DQ1 A(L» U=1 -Uf‘20 

03430 IF ITHASSILI .hEl*D«J0) NW=NVV + 1 

03440 DD 710 KAPG= I »MOANO 

034S0 IF tAlL«ArtPu) .NE* O.UDo} NN2TKH=NN2rR4F 1 

03460 710 CuNT INUE 

03470 C 

03400 C FClP »TEP~py-5TEP ANAuYSlS ADD THt HASS CONTft 10UTlDN-;i TO 

03490 C THE EUUATION DIAGONAL CUEFFICIENTS 

03S0O C 

03510 lFtNJVN»NE.4) GO 70 730 

. 03620 DO 720 l=l.NEaO 

03S30 720 All. 11 = AM»1) + A0« TMASSll} 

U3S40 kRin: 141 UAU* J}*l=l«NEad) iJs1,MSANl>) 

03560 GU TO 740 

03560 730 MR1TCL4) UAl l>3}<l=ltNEau)fJ = l«HOAND)*UdlI«l.) 1 l=^l*NE.Da)«L.=liUt.l 

03570 740 MRITE(9) t TMASSl I 1 « 1=1 t NEODl 

035BU C 

03590 IF 1 ,N0r. GENPHT) Gu TU 770 

J C3600 URITE 1u*200d1 H 

03610 DO 750 1=1«NLQ[J 

03620 IF IGENPCnt WRlTb 1 1 I , 200 7 J I A( 1 *3 1 . 3= 1 * HdANJ 1 

03630 750 14R1TC iu.AOO^U A{ 1 • 31 tJ=l tMSANDl 

j 03640 ttRITL Ior2009) 

Q36S0 DD 7uO l=l*NLQj 

03660 IF IGENPLH) MR [I £ 1 11 « 2007113(1 *3) 

03670 760 MRITE 16*200B] (»( 1« 3) *3=1«LL1 

J 036BO teRlTE (u*2010) 

03690 IF luLNPLH] hKIIE 1 1 1 ■ 2007 J I THaSS 1 1 ) « 1=1 «iScuu) 

03700 «A1TL lo«20Ucl) tTHASSt n « l=i »NE3dI 

03710 770 IFlH.dO.NdLUCK.) CO TO 1000 

) 03720 DU 700 L=KtNE2b 

03730 ANUUrt=ANURM + A(L*i) 

0 3 740 IF 1A(1..11 .NL. OtUUO .aND. A(u* 1I •LT. ANINJ ANlN=At(.» 1 ) 

03750 IF 1A(C»11 .Gt. hMAAI AHAX=A1L>1J 

) 03760 IF t All.* 1 1 ./4C..O.U0J NDEG=NDEG F 1 

QJ770 IF IA1L*1 ),tU*0»D0J A {1- * 1 } = 1 »Dl-20 

03760 IF { THASsIL 1 .NL *0 .001 NVV=NVV + 1 

03790 DD 7bO KA=l,rtMANO 






r~i . rzn [ZZD dJ dd Ld 




*F-Cr* 


PURTRAN V.oAtbf'i) /KI ^J-rtAlf-7U 


PA(jc 1— J 


03BQ0 

oaaio 

03320 
038JD 
03840 
038S0 
03860 
03870 
03880 
03890 
03900 
03910 
C3920 
03930 
03940 
039S0 
03960 
03970 
03930 
□3990 
04000 
04010 
04020 
04030 
04040 
04050 
04060 
04070 
04080 
04 09 0 
04100 
04110 
04120 
04130 
04140 
04160 
04160 
04170 
04 180 
04190 
04200 
0*4210 
04220 
04230 
04240 
04250 
04260 
04270 
04280 
04290 
04300 
C4310 
04320 
04330 
04340 
D43&0 


IR lACLftiK) 6.000 J NR2THH=Nr«ZTRfH-l 

780 CDNIINUE 

lf^lNDYN.NE.4l 00 TO 800 

DO 790 I=A«r4828 i 

790 At 111 I ^ An*!) + A04 TMAS51X1 

tiRlTEC4) UAfIiJl*l=i<tN828}.3=l (rtbHNOl 
GO TO 610 

800 teRlT£l4J ( (All >31 * l’*:K)M£2D) tJ = lf H3AN31 li ( 8£ liLI t 1=A»Ne:.28) »E-1 tl.L.) 
810 IfRlTE 19} [TMAS5111 •i=AiN£28) 

IF £ .not. GGNPRIl uU TO 840 
HPi=H+l 

KKITE (6i200o) MRl 
00 a<0 l=AtNE23 

IF IGENPOHl WRITE Ul it2007H All ■ Ji >3=1 tHOANUJ 
820 WRITE 16,2008) (All* J1 tJ=lfH3AM0l 
WRITE £6,2009) 

DO 830 l3K,HE2a 

IP (GENPCH) WRITE ( 1 1 ■ 2007) Onl i >3| ,3=1 ,1.EJ 
840 WRITE £6,20081 (BCl, 3J ,3=1, LL) ' 

WRITE £6,2010) 

IF IGtNPOrt) WRITE £ll,2007)(TMAa6(I) ,I=A ,.Nc,i8J 
WRITE £ 6 , 200 a)( THA5S( n , 1=14,NE2 b) 

' 840 IF £MH.EQ,HB) HM=0 

lOQO N5hiFT=NSHIFT+NE28 

IF (ROEG.GT.Q) GO TQ 850 
WRITE £6,1010) 

STOP 

850 ANiiRH=(ANQRH/NUEG)»l.U-8 
N r ERH=NE8 6*HQLaCl< *MaAN3 

PCT=iaO. 0D04DFU£JATtNN2TRM)2DFLJATlNTERMJ 
RATI 0=1. OD30 

IF (AMIN .NE. 0.000) RAT IO=AHAA/AH IN 
AAVG=ANURH41 .008 

WRITE £6,2020) AH iN, AKA2 ,HAT 10 , AAVu>P*-T 
RETURN 

1002 FORMAT £4FIO.O) 

1004 FORMAT £81S,3FI0.D) 

1010 FORMAT lolMOSTHUCTUKE WITH Nu JuGREES 0= FRkEOOM CHEEA OiiTA J 
2000 FORMAT (2/^ tOH STRUCT JRE ,13A,7HtUtMENI . 4X ,4HuUAO ,4X, 

A tlHHUl.TiPElERS,2 lOH t-OAD CASE • 12A, Im, 9X ,1 Ho ,9X , IHC ■ 9X.1H0.7 1 
2002 FORMAT ( 16 , ,4F 10 .3) 

2004 FORMAT 146HIS TEP-8Y-STEP :aOi-UTiON , 

X 37HC. ONTRQi- I NFORMATILIN, /// 

X SX, 35HNUH8ER UF IlMfc VARYING FJNcTiUNS =, 16 // 

X 5X, 35HmROUNO MOTION INOl CATOS =, la / 

X 8X, lOHEa.O, NONE, / 

X OX, 29HGT.U, RtAi) ACCEcERAT ION INPUT, // 

X 5X, 35HNUH8LR OF AHRlVAC TIMES =• lO / 

X 6X, 26HEU.0* ALC FUNCTIONS AHRIWE, / 

X 8X, 18H AT TIME 2ERU, // 

X bX, JSHNUMtlEK OF SOLUTION TIME aTEPS =. 

X 5X, JbHUUTPUT £ PRINT) INTERVAL =* 15 // 

















ApaSTF LAbH^.FUh 


f-PHTHAl* V*&AtOZlt /<l 


lai^a i-4 


0A360 
P 4-3 70 
04330 
04390 
04400 
044iO 
044a0 
04430 
04440 
04450 
044^0 
04470 
04430 
04490 
04500 


■4 aX » 

4 5/. * 

4 

A 5Xii 
« 5X» 
2D0a FOFiHAT 
dCOit POKHAT 
Jt 

*007 FORMAT 
2Q0ts FukHAJ 
2jQ9 FORMAT 
3010 FORMAT 
3030 FORMAT 
X 
X 


-* ul4#4 // 


=» bi4«4 // 


= c iil4.4 /// 1X1 


35«»0LUTT0M TIME INt-RtHEhT 
3UXHA5:;- FROFORTlOrMt, OAHFlNiS* 

JbRCoEFFlClttJT lALFi^Aj 
JOHSTlFFROSS-PROPOfmOWAO OAHPiRu» 
J5riLOeFFXCl£«7 tBETAJ 
tI6«7X»4Fia.^> 

tUOVERAOO MATKiCES* uLJC4« > 139/3 ■ 

• STIFFNESS MATRIX • J 
ttiP8E10»3J J 
i/t £1H 91PI0013.4T) 

T///j* LtfAO MATRIX *J 
K///t * MASS MAT*ilX *» 
t'OSriFFMESS MATRIX PARAM&rcrfG»»//9 
lPX,<MlTflMUH NOh"Zi;BO OlAOURAu cOEMENls J , iPoiO- ^t/. 
lOXt'MAXlMUM OlAOaMAO EOCMERT =•. OIO.S*/. 


r- 



04S1O 

X 

1 OX.^MAXIMOM/MINIMLTM 


= •* 

DlO.3 

»/■ 


f 

04S20 

X 

icx ,» average 

OlAGURAU ELEMENT 

-*# 

010.4 

f /* 



04S30 

X 

16X.*OfcrtSlTy 

UP THE MATRIX 


6pF5 . 1 . 

» PCT.+J 



04S40 

3000 

FORMAT , T3 7R0444 tfiRoR ZERO T tMt 

STEP. / 

1X1 



i 

04SS0 


ero 









COMMON OLUCXS 









f . : 

/DyM/t+14J 










»T 

♦O 

NOT 

+1 

ALFA 

4-3 

OT 

+4 

BETA 

+6 


NF« 

-no 

NOW 

♦1 J 

NAT 

+13 

NOVrt 

Flo 



1 

/EM/ 1 HU 70 i 










LHD 

FO 

RJ 

■H 

EH 

+3 

IPAO 

+101 

SS 

+103 


/EXTR/I/1+3J 









1 

MODEX 

+0 

MB 

+1 








✓crL/t+ri 










PRAV 

♦ 0 

EuPHT 

♦3 

ELPCH 

+3 

SbNPRT 

+4 

GENPCit 

+5 

i 

DEFPtH 

+4> 










SOOPMOOHAMS J-AouEO 








J. ; / 

OFI-OT, 

3S0RI* 










scauars 

AMO AMRAVS £ F4*f RU 

EXPElCl f ilLF Iti 17 1 4fi 

- tiut RcFEMEMEcP 4 


) 

«S0030 

1 

♦TETA 

2 

.10011 

4 

♦6«_ 

4. 

♦amin 

6 


«IP010 

10 

♦ 33 

tl 

♦NEQLi 

1* 

♦ 1H5 

14 

♦071 

14 


4HTEBH 

to 

.50047 

17 

•S0U4S 

30 

9 SPO 4s 

31 

♦NO60CK 

3* 


•50D34 

23 

♦II 

24 

♦NOES 

35 

3TR 

36 

•SU044 

27 

t 

950033 

30 

♦A 

41 

.503 J 1 

43 

. 9oO^^O 

44 

0 

34 


♦RATIO 

3S 

•S0P47 

37 

*NA 

40 

4AMAa 

41 

.50046 

43 


♦RUHEL 

44 

« 5U04*» 

4» 

♦nnztrn 

46 

9S0J4>« 

47 

.50044 

bu 


.50042 

G1 

4t>ri? 

S3 

•50041 

S4 

9 s 0043 

»S 

♦ANOMH 

S6 

i 

♦ RE20 

OP 

♦WijMJpr 

SI 

♦MB 

63 

.SOOs / 

U4 

♦ M 

64 


•SPPS6 

6y 

•SOOSs 

S6 

♦ j 

67 

♦MM 

70 

•50064 

71 


*SCPS3 

73 

.50007 

74 

.S0')j3 

74 

♦AX 

7s 

.50006 

76 












IZZ3 LIZl ■ CZZ] czn CZ3 c 




cz:} c 


ADDSTP 

EASHEU. 

FOR 

FORT RAM 

V*oA( t^El 

1 y^i 

24 -rtAr- 

75 

i 06 

PAOE 1“6 

.SOUbl 

77 

.60005 

100 

•sooso 

lux 

• 50004 

X 02 

♦ 11 

103 

>30003 

104 

A 

105 

■50002 

106 

•50001 

107 

♦NVV 

ilO 

.soooo 

111 

♦i-HN 

112 

«NUMb 

113 

• 50067 

U 4 

•S 0066 

1 15 

*50066 

116 

.S 0064 

117 

♦HPl 

120 

•S 0063 

121 

* 500 X 7 

122 


133 

•50064 

125 

•500 16 

126 

per 

127 

.50061 

140 

•soois 

131 

*AD 

142 

•SOObO 

134 

4 KAP 3 

145 

.50014 

146 

. I 0 Q 07 

14 7 

•50014 

140 

•10006 

141 

4 NUH 7 

142 

•50012 

144 

•10005 

144 

• 50011 

145 

•10004 

14 b 

•SOQIO 

147 

•10003 

150 

* 1 . 

151 

. I 0 DO 4 

152 

4 MUAMU 

lb J 

.10001 

154 

♦AAV 6 

156 


157 

•IQOOO 

160 

TMAS 5 

161 

4 Nc 55 

162 

•S 0 Q 27 

164 

.50036 

164 

.50071 

165 

•50025 

io 6 

4 MSHF T 

167 

.50070 

170 

.SQQE 4 

171 

•50043 

172 

•S 0022 

173 

•50021 

174 




TEMPORARIES 
.AOOIG 616 

ADOSTF X NO ERRORS UETECTED I 


O Q 

gi 


V — J 

03 


£ 

I 





r Q4±>60 SUtiUOUTfNd (.AUzlANmaANJiNalFjLHtXMiSfr^tNtlaNUMtNbJ 

0Ab70 IMPLICIT *^fcAL»a t A-rt, U- jI ) 

O^SHO LUbILAL uLPRT. cLPCrlf btiNPRT • btNPCrt» JtFPCri 

04590 C 

) ' 04600 C CALL bP BY RUaSi I£AH«PLNAA»BRlCAti»TPLAr£f CLAMP.CLSTaO.PXPEK 

04610 C 

04620 C- CALCULATcS oANJ tSlJTri AND kfRITCS 5TIFFM565 HATRIK UN TAPE 2 

04630 OIKENSION LM(U «XM(U NCiHTaPiNDHtA) 

' 04640 COMMON /EXTRA/ MODEXuNTB 

C4650 COMMON /CTL/ GRAV»£lPRT »ELPCH.«ENP ttT*OcMPCrtf OCFPCHt 

04660 1 iMTCGtMlNJNO, IPLT 

04670 HFi=2 

' C46flO IF (MlNBNO ,ufc. 11 ^Pl=7 

04690 NF2=1 

047QO IF (HINOWO aGE. IJ NF<*=7 

04710 MIN=ilOOOOO 

> 04/20 HAX=0 

04730 00 500 L~1«ND 

04740 IF fLMILlaEQaOJ GO TO dOO / 

t 04750 IF li.MtL}.GT«HAXl HAX=LMIL1 

i ( 04760 IF (LH(L1 aLT.MlNl MlN=LMtt-l 

04770 aco CUNTINUe , 

i 04760 NOIF^MAX-MIN+1 ■ 

04790 IF INDIF .LTo 01 N01F=0 

! J 04800 IF INUir.GTaMdANDl MbANO-MOlF 

h 04510 IFtMOilEXaEQalJ GO TU 51u 

i 04820 C ■ 

i 04630 LRD=NO+tNO+ll/2+54ND 

I 1 04640 WRITEINFI) LRO t NO ■ (LH( 11 « 1 = 1 <N0 3,f ( (Si I > J 1 « J=i tNUl • 1 aNO 1 • 

f 04850 X t(P{l>Jl il^ltMOl .J=l«4l9(XHm aI=l«N01 

04060 IF («NUTa Ei-PKTl RETJKN 

■ 04870 WRITE (6.1000Mc-HUlal=l«N0) 

I 04860 WRITE (6.1010) 

1 04890 DO 900 1-l.NO 

[ 04900 IF (ELPLM) WHITE U 1 1 1020 H5 ( i t J1 a 0=1 »N^ 1 

04910 OOO WRITE iOiIOSOHSl IiO) .0=i«ND) 

; ) 04920 ■ IF (ELPCH) wRtTE U 1 . iOalO J { i PI 1 a J 1 a 0=1 ■ 4 ) a 1=1 • NO ) 

C49J0 WRITE (o«1040H(P(l«J}a0=l94}»X = l*N01 

04940 IP (cLPCH) WRITE U 1 « 102U )( XH( 1 ) a 1=1 1 NO ) 

04950 WRITE ( 6 a 1 050) ( XH( I) a 1=1 «ND) 

J 04960 RGTUHN 

C4970 aiO WKlTEiNFE) >Nl>« NS a ( LM{ X 1 a 1 =1 a NO) 

04960 RETURN 

04990 C ■ 

) 05000 1000 format ('OELcMEMT i-UCATION matrix °a/a(iH alUiio)) 

06010 1010 FORMAT (•OElEMENT STIFFNESS MATRIX ■) 

05020 1020 FLHHAT UiP6EI0a31) 

05030 1030 FURMAT (/a(lii alPiOcl3-4ji 

J 05040 1040 FORMAT 1'OELEMENT LOAQ MATRIX ‘a/ftlH alP4uI4a4J) 

0SD50 tUbO FURMAl ('OELEMENT MASS MATRIX 'a/a(lH alPlOElJaA)) 

OsObO EMU 

) ■ 

COMMON QLUCKo 


/ExTRA/(fr2) 




CALfJAN LASHEi. ».FOR 

Ncocx +0 nr«i 


FORTRAN \Z-5Mo;-lJ /Ml 
+l 


XJ-MrtTf— 7ii 


iz;j8 


f-»Auc 1-1 


/CTU/'t+J21 

uRAV 

OEFPCH ^6 


fcuPST 

XWTCG 


■i-ti 

♦7 


ELPCrt 

MlNiJMO +10 


OfzrtPRT 
IPut +11 


(jEMPCH +£> 


SUBPRQGRAH5 €Al.t.£D 


SCAUANS AND ARBAVa 1 *•♦« rtO EXPLICIT OEplNlTlUN 


«UT BEFtBcNC&a i 


C; 

< 


.50020 

t 

»HIN 2 

4NS 

J 

4NF 1 

4 

♦MAX 

b 

*ND» 

6 

4ND£F 7 

♦NF2 

10 

5 

n 

P 

12 


13 

.50007 14 

• 50(106 

15 

.50005 

16 

.50004 

17 

.SOQ03 

20 

.50002 21 

.50001 

22 

.50030 

2>> 

•sooir 

24 

•S0016 

2S 

.50015 2o 

Art 

27 

•50J14 

30 

4rtD 

3l 

.50013 

32 

.50012 33 

.£0005 

5* 

•soon 

35 

•10004 

3d 

■5C01O 

37 

.10003 40 

4l_ 

41 

•10002 

42 

4M5AN0 

43 

.10001 

44 

#1 45 

•lOQOO 

46 

UH 

47 

.50021 

50 

91-RO 

bl 








T£HPCRA(UE5 








.A001& 

122 








CACBAN 

1 NO 

tBROBS JLTcCTHO i 








1 

\ 


I- 


T 


LZ3 


CZJ : illlj 


LZ2] CZl] □ C 


C™] €Z3 1Z13 LZD dH d3 dI2 [1123 C 


|_i-A 


CHECK LASHcL.FHR FLiftT«AN V .SAt £>2i ) /K1 £:j-«AV-7d l^i^o PAut i 


\ ' 050 ro 

i C506U 

I 05090 

I 05100 

t t osno 

05120 

I 05130 

I 05140 

) 05150 

t 05160 

05170 
osiao 

: • 05190 

; 05200 

> 05210 

i 05220 

i 05230 

1 . 05240 

I CS230 

; 05260 

05270 
05250 
05290 
05300 
I 05310 

05320 
C5330 
05340 
' 053SQ 

05360 
f 05370 

t 053S0 

; J 05390 

0S40Q 
05410 

i , 05420 

i 05430 

L 05440 

05450 
OS460 
1 05470 

054U0 
05490 
05500 
) 05510 

05520 
05530 
05540 
) 05550 

OSSbO 
05570 
05500 
> 05590 

05600 
05610 
05620 

i , 

V 

•oj 


SOHRC J T 1 ME CHEC K n D 1 1 FwG t NEu t N J MciX iN EuT r H 1 
IMPLICIT HEAL45tA-H#U-ZJ 

LCOICAL ELPRT.EL-PCritviENPKTtbehPCrtfOEFPwri 
C 

C CALLEO 3Y MAIN 

C 

C THlb dU3Kt)U7 1NE CHELKS TO SEE THaI Auu 3Euri£Eb OF FREEDUH HAVE 

C STIFFNES5 ATTACHEU 

C ■ ■ ' • ' 

CUMMCrt /ELPAR7 NPARU4J uNdMM? 

DIMdNSlUM lFuC{NEO>»ldtNUHNP*ol 

COHMOH /'EM/ LRDvNOiLMCbdJ t lPAOt534c^dil 

COHHUM /CTL/ ORAVacLPRT f fcLPCH.yEM.’RT . ^sNPwH, UcFPcrt* 

1 lUTCvJsMlhLNDf IPLT 
COMMON /EXTRA/ MOOEA 
NF=2 

XF (HJOEX acU* IJ NF-1 
Q 

L SET THE CHECK VECTOR TO 1 

C 

DU 100 1=1 .NEQ 
IFL«tU=l 
100 CUNTINUE 
C 

C READ THE CM VECTOR FUR EACH ELEM&NT AND 5ET THE CHECK VECTOR 

.C 

c data check mode is off 

c 

IF iMiiDEX >£Q« n ^O TO 200 
REWIND NF 
DU 120 N=ltNUMEL 

READ INFJ LROtND«(uMl 11 a 1-1 tNJJ *15511 1 *1=1 aLRiJ) 

DU 1 10 I=1*ND 
IEQ=CMm 

XF (lEO *CE. 01 GO TO 110 
lFLutlEQ|=0 
110 CONTINUE 
120 CONTINUE 
REWIND NF 
GO TO 400 
c 

C JAIA Cni^CK MUdt: ON* NU dAMJNiurrl MiNiMl/AllDN 

C ■ ' . 

2UQ IF IMINUNJ .GE« 11 GO Tu 30U 
REMIND NF 
DO 230 N=1 *NELTVP 
READ INF! NPAR 
NUHE=NPARt2J 
DO 220 J=1*NUME 

READ INFI NDiNSslLNl n *1=1 <Nl>1 

DU 210 X=1*ND 

IEQ=LRt[} 

IF ( IeQ .lE. 0) GO TU 210 
1FLG{1EU)=0 
210 CUNTINUE 
220 CUNT INUE 






CriECK uAa):tEu«PO({ 


FJfURAN V^5A(6<>i} /M -7tJ 


X- X 


05630 

260 

caHTlilUE 

asfifto 


REWIND NF 

C56aC 


iaO rO 400 

05660 

C 


05670 

,c 

DATA t-HfcCK viaoc ON «lTrt aASDwlJTrt ,<tl Nl Ml 2/tT (UN 

05660 

■ C 


□56S0 

300 

REKUND NF 

0570C 


UU 330 N=1.NHUTVP 

0S710 


READ INF) NPAH 

05720 


NUME=NPA»U2) 

0a730 


DO 320 J=1«NUME 

05744 


READ (NF) NJ«NSp£l.M(l} •I = 1«N31 

05760 


DU 310 l=ltNO 

C576C 


IEO=UM(U 

05770 


IF (lEQ .4E* Oi uU TU 310 

0575C 


IFLDUEQJ=0 

05750 

310 

CONT INUE 

05 BOO 

320 

CQNT INUE 

05 a IC 

330 

CUNTINUE 

05620 


RctolND NF 

05630 

iZ 


05640 

c 

SCAN THHOOGrt THE CHECP. VEt-TOR 

05650 

c 


OSBoO 

•♦70 

N=D 

05B70 


DO 430 1-iiMEU 

OS 660 


IF *£Q. 0) 3D TD 430 

CSB5C 


N=1 

05900 


HDDEX=1 

05910 


oa 420 N=1»NDHNP 

05920 


DD 410 J=lib 

0595 0 


IFtlOCNpJl.NE.l )6U TD 410 

Q5940 


t(RITLt6,2C0ajl ,N* J 

05950 


DU TU 430 

059o0 

410 

CUNTINUE 

0o37a 

420 

CQNTiNUe 

059B0 

430 

CUNT INDE 

05990 


IF (N .EQ. 01 WRITE IbpEOlO) 

C600G 


RETURN 

06010 

w 


06020 

2000 

FURMAT( • OCQUATIUN* ,15,' tNODE^iljt* 0 .F * • > 12 . * ) 

06030 

1 

*• STIFFNESS. CDnTINDH 1N THE DATA CH23jv MUDc*) 

06040 

20 10 

FDKHAT t*uALL CQUATIUNd HAVE SflFFNcSS AlTAcrtaD*) 

OoCbC 


END 

wCKMuN 

Bt-vlClvS 


/EUPAK/ 

t + 1 7) 


NPAR 

40 

NDHNP +16 


HAS NU* 


/£M^1 +l li 70 J 
1_RD +0 

uRAV t-J 


ND 

+1 

EM 

♦•2 

IPAD 

+ 101 


+ 102 

eePRT 

iwrca 

+2 

+7 

El-Pit,rt 

rUNDND 

+ 3 
+ 10 

muNPHT 

IPLT 

+4 

+11 


♦ 5 


©ft) 

G ^ 

'^05 


c;:r] L 




'CZ2- 


f ■ 



CHECK UASHE-.^Uk FoKTf<AN V.:jA(o^XJ /M <ii-,*4AY-7a liil^li P«v*c 1-ri L 


/£XTRA/t+l i 
RQOEX 4-0 


SUaPRUCRAMS C«.LL£k< 



SCALAHS ANJ ARrtA*» j 4. »'*'« NU EXPClCiT tlEFIMlTlHN - " V» .NJT Kc.Fi.fi fcN^(-.J i. 


4 = NS 

1 


2 

ANCQ 

3 

»N 0,4=.. 


♦WF 

5 

JFJ.G 

d 

lo 

7 

4^3 

10 

•50007 

11 - 

•SCCCu 

1= 

• SOtlOS 

13 

• 3QuJ*4 

lA 

•50003 

1 3 

• 50 JJ 2 

is 

•SuOOl 

17 

«S0000 

AO 

3NUMt 

21 

.50017 

Ad 

*Ncl-T Yp 

23 

•suois 

2A 

•scois 

2s 

-SOOl^ 


•50013 

27 

. .SOOl^ 

3C 

•soon 

3i 

e JOOOA 

32 

« 50010 

33 

•10003 

3A 

• ii:/ji]2 

35 

•10001 

3& 


37 

■ 1 0 00 0 

AO 

*1EQ 

41 





TEMPCRAKIES 









.AOOlO 

100 










CHECK JL JETtCTtia J. 



1 

E 





r" 1 . 


□ CZj 


Eizii czD czn izzi cm cm 


fcL TVP£ LA5 h£«_ o F Ok 


FURFRAN V,5A(0<ii> /Kl 




00060 
000 70 

06000 i, 

06090 C 
06100 C 
COHO C 
OOIOO 
06J30 C 
06140 
06150 
06I6C C 
06170 C 
06150 C 
C6190 


SUL.KUiJTliF£ (_LTVH£lrtTyPe> 
lMFl.lCir HEAL*6(A-H«a-ZJ 


CALLcO 67 MAIN* STRESS 


IF KMTVRE-IS) »ca. 0) 6U T O 13 

IF {imTPB-i4.y .SQ, py (sp to /3 

11 pfrtlTi- 16*100 1 MTYPc 
G6 7L 900 


06160 C MUl-TIPtd i-AVERcTj tlYI-lKOERlCAt. SHct.1- 

C6190 13 Call t-AYShL 

06200 900 RETJKN 

06210 C 

06220 lOJ FUH.MAT t 'OE lEHENT* » 14* • IS NOT IHPLEMdNTED VET*) 

06230 End 


suoPHgsrams called 


LAYSHL 


SCALARS AND A4BAYS i "♦« NU EXPLICIT DEFINITION - NOT K_FERENC£0 1 


♦MTV PE 1 


TEMPORARIES 


■AOOlo 13 


ELTVPE X EHKCRu uETcCTED ± 






•Mm 




•msm 




2J2|fg2j2 




ERROR 


i.A3Htx*FIjH 


fORTRAM Uk6A(b2i} /Ki 




Id: PAbc 1 


062AU SObRUJTlNE eRRijri(R) 

06250 CCiMRQR /cLPAR/ WPARt 14) » NJrtWP.M6AN0*Nt£Lr VP*M,Nt;. NJiM^^NSiiRTtJr^NL- J 

t)626C) IF tN *GT>, MTQT) GO TO 100 

06270 1»RITE I6*2000) N 

062Q0 RfcTURN 

0629C 100 liRlTE t6i201Ql N«RTOT 

0630 C STUP * ' 

06310 2000 FORMAT (• 0RL.QUI REO BLANK LUMHQN FOR TrUS STbP=**17J 

06320 2010 FOKrtAI ( '□♦♦♦FATAL ERnOK#*** «/. 

06330 A* REGUIREO blank COMMON =*, 1 1 G* / p 

€6340 X* AVAILABLE BLANK CUMM 3N=* , 1 IC ) 

06350 END 


COMMON BLUCKB 

/ELPAR/t+3in 
NPAR FD 

N2 +22 

NEQ +27 


NUHNP 

N3 


MdAND +17 
N+ . +24 


NcuTV? +20 
N5 +2.5 


SdUpRCGRAMS CaLLEO 


5CALAH6 and arrays i NO EXPLICIT DEFINITION - Av4'» NOT RbFERENLED i. 


TEHPCRARICS 
aAO.Olo 3a 


ERROR 1 NO ERRORS BETECTED 1- 


f 


‘s 


r 




CHj C 


CZ3 !~J 




■E^ B X 


HEDPRT i_A3Hcl.*FuH FUHTKAN V*a,«iCti21) /KI *iJ-Mn»-7o 


PAOE 1 


‘ . 06J60 

C 6370 
U63a0 
06300 
06400 C 
06410 C 

C64a0 c 

06430 
06440 
06450 
06460 
Co 47 C 
06430 
06490 
06500 
OuSlO 
I 06530 

06530 
06540 
0655 C 
06S60 
06570 
C 656 C 
06590 C 
\ 0660 C C 

06610 i. 
06620 
C 6630 C 
I 06640 

06650 
06660 
06670 
I 0663 C 

06690 
06700 
0671 C 
t 0 a 720 

06730 
G 6740 
06750 
p 06760 

06 770 
06730 
06790 

} 06300 i. 

06 CU 0 C 
C 6320 
06330 
I 06340 

06 350 
0686 0 
06870 
) 06830 

06590 
0 o 90 C 
06910 

1 


SUBrtCjriNu rttOPKT t( 4 urN) 

IHPi-lClT 

REAU 44 0 AYTMpD 3 ATE 
REAU *4 Tf St SbRAVtHdJ 

called by main 
COMMON /uAR/ A( 9999 ) 

LlJbl CAL ElPRT t ELPLrtt ocNHKT p 66 .NPCii pDcFPwrl p 3 EOS I 
LUHMON /junk; /HbU( l 6 j tOJK.( 400 ) 

COHMCN /^.lPAR/ NpARtl 4 ) p NUMNPpHBAN 6 «Nct.r rp ph 1 p N 2 p .>( j pN 4 pN 6 p MT UT pNcQ 
COMHJN /CA / uaUlt:bbJJ 

73 YN 7 IDuatll J pNJYHpT^□PLr|pI^i 3^_^1 
COMMON /TAPES/ NQQlb) 

COMMON /EXThA/ MOOEX p NTb«N 10 SV p NT 10 t aEUB pNUMEL p T I 12 1 >S 1 12 i 
COHHuN /SOL/ NBLDCXpNEaHpLCpNFi 10 Lm«NElSpNADpNVVf ANORHfNFO 
C OMM DN / CTl/ GK a V p ELP RT • ElPCHp GcNPHT # o£h PC ri » UEF PCH p I WT Ctt p 
I HINBNDpIPlT pNPLUTSpNPCT 

COMMON /GEOsTF/ GEOST»NeUGEO 
COMMON/'SElS/NRSCpXAPGpNRSp IlDUMpODATElE) 
cummon/fokhat/na 
lUMMON /JAYTIH/ OAYTMtS) 

OlHtNSia.s IPRT 15 ) 

PROGRAM lAPACI TV LGNTRQlCEJ 3Y TNb FOLi^JmING TWO STATEHEnTS 
MTOT- 9996 


6 ftdVitNO I 
RElklND 2 
KEMlNO 3 
RLttlND 4 
RCMNU 7 
KtSlNO 9 
RLmINO 12 
HEttlNO 13 
RLlmiNO 14 
NTb - a 
KlsINO NTb 
NTIO- 10 
rtSWlNU NT 10 
NELaEO=0 
Nl=l 


P K U o H A M 


L □ N T R a L 


3 A, T A 


head tap lOOpts«J=99C ) rttiO p NuMNP, NECTYP n_t. pNF p NJYN t MDD£X,NAJJp 
I KEObpMlOSO p IPRTt 1 p(TC« p MINGNU p IPLT pNRSCpbRAV 

NA-NUMNP 

IF lnumnppEGpO) go ru 090 

IF tuRAV *EOp Op JO} CKAV=360,4JU 

1 
Z 
1 

2 
1 


eCPRT -IPHTlal *GE. 
ELPi-H =1PRT13) pEO. 
6UNPRT=lPKTt4) .Gbo 
6ENPLH= I PRTt 4 ) *Eu» 
0EFPLH=1PRT lA J tEQ. 





3 '**1 










r^TT^r 


tAb At 







rlEOPBl 


Qt9'ji0 
06930 
C 6 ^ 4 & 
069S0 
06960 
0b970 
C696Q 
06990 
07000 
CTO 10 
C70£l0 
07030 
G 7040 
07050 
C 7060 
07070 
GTOaO 
07050 
07100 
07110 
07120 
0 7130 
CT14C 
07150 
07160 
07170 
07130 
07150 
07200 
07210 
C722C 
07230 
07240 
C7250 
C7260 
07270 
07260 
07250 
073CC 
07310 
07320 
0 7330 
0734 0 
07350 
07360 
07370 
073SC 
07390 
0 7400 
0 7410 
07420 
07430 
0 744-0 
07450 
074ei0 
07470 


i-ASHE^*FC/l FJRT)l4ri V.SAloSlj /Al 2J-H47-76 12i3o PAiiE 


IF (HUJ6X«6T«Ol rtOOcA i 

C ' 

C bill 7i. 0L,uCKt-i£TT6« HkAUlitG FUR THE «E* H203ucH 

C 

C wAlTE 16,20001 

C VHITS. IStStOOU 

L WRITE {6,3S0J 

C 

C OiET THE CURREMT DATE AMJ T 1H& 

C , 

CALL DATE 
60RAV=SN2l t uRA Vl 

WHITE 16,200] rJcO,HtrHHF»HELTyp,Li-»Hf' 

WRITE {6^2011 /lOO£X«HAU,K&Qb^MlU6V,5i,RAV 

WHITE {6,210> MTOT 

1F(XE03,lT.2} K£QS = 99999 

IF INOyM.NE*OJ LE=1 

ZF£i-L*GEilJ Ga TO 10 

WRITE 16,3001 

STOP 

C444 DATA HURTrtUi.e 3AVE 
10 IF (HCOEX*cQ*li 

AWftITE I/JT3J ricD,Nu.<HH»nELTyFft3.,?JF,ftLjyr( 

IF IH1N3H0 .GT# OJ WRITE 16,370) 

IF tiPET -GT, a) WRITE (6,380) 

1F1XP1-T-LT,0 JV(rtXrEl£.,390a 
C 

KJVrt - lASElHOyH) 41 
CEOaT^.FALSE, 

IF (KDYH.LE,6) HETtJhM 

nuvri=KDyr(-io 

Gi.05T^.TRlJE* 

WRITE 16,360) 

IF IKpyK ,LE, 5J HETOKN 
WHITE (6,310) Mjyf# 

990 4.U)4T1I4UE 

1 F ( «HL JTa , uT , 0 ) STOP 
REWIND 3 

WRITE 131 LL,HF ,AuyN*H60d,N£U,Nou6LX.,WJjE4,rt0HNP 
. STOP 

too FORMAT (l6A4,XfOl£i,6TI *^l5,la,Hl0*0) 

200 FORMAT (IHO . 304* 13A4, / . 


A 

36mc O « T R 0 O 1 N F O R 

H 

A T 


1 0 ti. 

X 

27H 

NUHdEK OF HOOaL POINTS 


1^ 

/ 

4X, 

X 

27M 

HOrtuEK OF element lyPES 

s f 

16 

✓ 

4X, 

X 

2 7H 

NUMdEH OF LOau LASE5 

2T f 

IS 

/ 

4X, 

X 

27H 

NtlHiiEH OF FREOOENlIES 


15 

/ 

4X-* 

X 

27H 

rtNALr:>IS coot iNdYM) 


la 

/ 

4X, 

X 

lorl 

t.d.0 , wTaTIl, 



✓ 

4X, 

X 

26H 

cO*i* MODAL ExrHAi-TlUN 

a 


✓ 

4 X, 

X 

25H 

cu,2, rUHCEO HLSPONSE, 



/ 

4X, 

X 

27H 

Eo«j* response spelthom. 


✓ 

4X, 

X 

2db 

EU.4, DIRECT INTESKATION, 


/ 

44, 

X 

4dH 

ED«s> FREOUENCy RbSPUNSc.) 




201 i 

F-ORHAI 

I4X, 





X 

L/H 

SuLOTICN AUOE IMdOLX) 

^9 , 

la 

/• 

4X, 



rteaP(ir uashel.for 


FURTHAN V.SAloiin /kI iij-(4AT-/<j 


111* ^6 Prtuii i-i: 


CT^iC 

07490 

07500 

07510 

07320 

C7530 

07540 

07550 

0756C 

orsro 

0756 C 

0759D 

07530 

C751C 

07620 

C7630 

07640 

074SC 

C7660 

07670 

07560 

C 7690 

07700 

07710 

07720 

07730 

07740 

07750 

07760 


^ lUH hU*Ct EKECUriUiH* / 4Xi 

X 2011 bO*lt uATa.ch5i.Ki / 4X« 

X 19H HUHlSlsR UF SUd6PAC£* / 4X, 

X 27H ITEHAriQN VECT0R3 CWaUJ =t 15 / 4X* 

X 27H eOUATi QN5 per BLOCK ~t 15 / 4X» 

X 27H TAPElO SAVE FcAO (NICSV) -i 15 / 4Xi 
X 27H GKAVlTATlGNAU CONETAMT =, ul2«5} 

210 FURHAT (5X.* TOTAL ULANK COMMON t MtOT , 15 1 

300 FORMAT t// 4dH *♦ ERROR, lAT _Ei»ST UnE cuAU ^A:jE 13 REuJlREOl 
510 FORMAT I// J3H 4+ ERROR, ANAcySlS COOe [NOYN =« Ioi 9H) IS 4A3. 
JbO FORMAT (//,45X»* USC VtRaluN JOulT 4# 1^7u* J 
360 FORMAT t»0GE0M£rRlC SUFFNEaS JPTIOM 15 BEING 05E0* ) 
j70 format i » OOANOwluTri M iN i Ml E hT 1 JN l5 REuJESTcu* 1 
JdO FORMAT t'CSTHUCTORE Pl_JTTlN« 13 KeOJESTbO* 1 
390 FQRMAU'OMESTART OF STRUCTURE PLOrTlN j REaUtSTEO'l 
2000 FORMAT { * i* */# 


— * ul 2 « 3 } 


X15X* 

X 

X15X> 

X 

X13X» 

X 

XL5X, 

X 

}il 5 Xf 

X 

X1SX> 

X 


2001 FORMAT 1 


SSSSSSSSS5 

vv 

SSS3SSS3S5SS 

VV 

ss ss 

vv 
ss 

vv 

S5S 

vv V 


SSSSSSSSSr 

vv vv 


AAAAAAAAAA 
VV* r/, 

aAaaaaaaaaaa 

VV* .7, 

«A A A 

VV * ly, 

AA AA 

VV • ,/, 

AA AA 

V • 

AAAAAAAAAAAA 

V • J 


PPPPPPPPPPP *il7X, 


PPPPPPPPPPPP* »17X, 


PP* • 17X, 


PPPPPPPPPPPP* ,1 7X, 


07770 

X15X.* 

SSSS5SS5S 

AAAAAAAAAAAA 

PPPPPPPPPPP 

»il7X, 

C7730 

4 • 

VV vv 

♦ »/. 




07740 

X15X»* 

SoS 

ah 

AA 

PP 

• *17X, 

07000 

X • 

vv vv 

* »✓• 




0761 0 

xlaX»* 

55 

AA 

AA 

PH 

* f 1 rx » 

C7330 

X '• 

vvvv 





07030 

X15K, >S£i 

SS 

AA 

AA 

PP 

« »17X. 

1 07340 

X * 

vvvv 

’ ./• 




07850 

X l5X« * SS5SSSSS5BSS 

AA 

AA 

PP 

* »17X, 

C7360 

X 1 

vv 

* 




07870 

Xl5Xi* SS35SSSSSS 

AA 

AA 

PP 

' . 1 7X » 

) 07880 

X • 

vv 

' J 




C7890 

2020 FORMAT 

1//,45X»5A4 

J 





COMHON dCOCXS 


/uAR/ 11-470 JQJ 
A +0 


/ JUNX2 (1-642) 
HEO l-O 

JUK 

+ cJiJ 






■ t. 

/ELPAR/ t-r 301 
NPAR +0 

NOMNP 

f Vo 

MbANO 

417 

NtbTltP 

420 

Nl 

+21 

N2 122 

NEQ 42T 

NJ 


N4 

♦ 24 

Nb 

♦25 

HTOT 

+26 j 







1N4. 


07910 

07920 

07930 

07940 

07950 

07960 

07970 

07980 

07981 
07990 
08000 
OBOlO 
08020 
08030 
08040 
08050 
08060 
0HQ70 
08080 
08090 
08100 
QBllO 
08120 
08130 
08140 
08150 
08160 
08170 

oaiao 

oaiet 

08191 

08192 

08193 

08194 

08195 

08196 

08197 
08200 
08210 
08220 
08230 
08240 
OB250 
08260 
08270 
08280 
08290 
08300 
CBJlO 
□8320 
08330 
08340 
08380 
0836 0 
08370 
08380 


c::r3 (zz 


B ^ cm cm cm 


cm cm c; 


cm cm 


LAaHEL.FQrt 


FDRTKrtM V*3A1621) /M 


aj-MAic-ra 


12 ; 30 




SU8RDU7 1 ‘til iNl-l 10t6rrMA3S«N0MN3*NE0iJfNS3 iLUaNuuuCh;) 
IHPi-lCIT R8Ai_*8lA-H,Q-Z) 


CALuEiJ QY MAIN 

1HP3T NO Li A I- LO AOS tSlATIC oJLUriONJ UH NOJAu NASSuS IpVNAHlCl 


OIHENSIUN lOtN^MNPio) * 31N30«1-L) *TMA5SINS8) 
COMMON 7UNII-30/ FP<(999| .FF71999* fFPZl999> 
COHMON /JUNA/ M l 6) 

COMMON /EXTRA/ MUucX 
COHMO.^/CTu/oKAV 


HEAP the IJ ARRrtf FBJH f ICc 3 
hU« 11 80 TO 105 


IF IMOObX 
:^pEVtlN0 a 
VaHEAP 181 IJ 


ZcRO THE MATRICES 


00 100 N=1»NS3 
THA36tN)=0 .000 

00 100 K=liCL 
100 a(N«A)=O.OOS 


REAP THE CAROS ITERMlWATt WHEN THE NUOc NJH3CR“0) 


lOS 

110 


140 


11 s 


118 


120 

130 


skRlTc: (CiEOOOI 
READ (StlOODl N 
Cfsl 

Rtl)=FFXtN) 

H12)=FI-Y(NJ 
K13J-FFZIN1 
R(41=0.D0 
ft(5J=0.D0 
R l6j —0.00 

It* (N .t_6. 0 1 oj TO 
IFlt. .8T .0 IGU TJ 118 
DO 11b 1=1 »6 

IF iRI I 1 .t-T.O.ODOlKt t l='rt£ I 1/L.KAV 
IF £Rl2).E0«0.0J0)RU)=Hi 1 ) 

IF IR13) .60.0.000 )RlJl=Kl U . 
CONTINUE 

WRITE (u.2010) N.C.k 
IF tc .OT. Uu) 30 TO 
IF IN .vjT. NUHNP) OC 
IF iMCDEX .EO. 1} CU 
DJ ibC 1=1*6 
lEO=linN, 11 

IF ( l6Q .EE. 0) uU TP IJO 
IP lE «eE. 01 oU lU 120 
bl 1 eQ.E1=R( 11 
OU TJ 130 
THAS6I lEoJ-iU 11 
CUNTXNUE 



20 0 

TJ 210 

ru 110 








^•■ faLr.Aj ■AV-.a , 







1 ■ ■ , 

IjMl- LASHtJL.hiiR PJrtTRAN V.3A(ojiU /M i-3“MAy-7o pAuii 1-1 


! ‘ Dd390 CD TCI 1 10 

{ 0 B 400 C 

0B410 C URITt Trie UuAU AND MASb' MATRlLcS UN h J 

0B420 C 

i B8430 140 IF IHQJEX .EO. 1} HETUKN 

r OB440 REMIND 3 

c 084S0 OJ 150 N5=ltN5LUCK 

E DS 460 Nl=NEQa«(NU-iJ+l * 

I i 0 S 470 N£=NEQ 54 NtJ ' 

I 08480 «fUn£ ( 3 ) UBtNiL)sN-Hl*N2)>4-li>Et,)itT.HA5SlN)»N=Nl*N<il 

' 08490 150 CUNT INU£ 

i 0 B 500 HtTUFlK 

1 oesio C 

QS520 C ERRUH nANt)LlN5 

ii' 0 8530 C 

I 08540 BOlr HRITE (5*3000) L»LL 

[ \ 0855Q MUDcX=i . 

t OaS6Q Gi) TO 110 

I 03570 210 URlTu (6*3010) N.NUUNp 

08580 HOOcX=t 

? 08590 Gu TC 1)0 

08600 C - 

08510 lOOO FORMAT (13) 

08620 41000 FORMAT ( 'IN 0 O A t. I. O A J 3 (5 T A T 1 C) O. ft ** 

[ t 08630 X'H A 5 5 £ S (DYNAMIC)' >///« 

' , 08640 X3X, ' NODE* *3X> 'LOAD' ■ 

■f 08650 X2l9X**X-AXlS»,9Xi'V-’A)aS'*9X**Z-AXI5'>*/*' NuMoER* <3Xt 'CASE* * 

, V ? 08650 X3(10X.'F0RCE'Ji3(SX.>HDHENT*)*e<'.1X) 

‘ r 08570 2010 FORMAT 1 2 ( 3X *14 ) > IpbElS* 5 ( 

08680 8000 FORMAT ( ' 0444 FATA 1 . ERRDR 444 LOAD N JMbeft* • 1 IS • ' IS GREATER** 

OB 690 X'ThAN THE MAXIMUM* • *1 ID ,✓, < COMTlNJe 1 '• DATA uHECK HUDE') 

08700 3010 FuRHAT C 0 ***FATAC ERR 0 R =«'44 NDJt NJMouR' , I 10 ♦ ■ tS UREATtR',. 

> 08710 X' TMaN Trifc maximum nUMcIER UF NJDE5 1 • . 1 1 J * * ^uhT INJb IN DATA'* 

08720 X' thtilN MODE') 

00730 END 

t ' . 

M 

I CCMHUN 8CaCK.5 

i - ■ - 

/UNIL 00 /l 4 -lJb 52 ) 

» FFX 1-0 FFV +3715 FF2 +7634 



/UUNX/(+I4) 

R +0 

J /'tXTRA/t + l) 

MUUEX +0 
/CTLXl+2) 

GRAY +0 

I 

SUDPHUliRAMs i,AUCtD 


1 





ltd. FiiUlHAU V.3Aler21J /M ^3-HAy-l'o i-^ 


i t 

i: f 

i 

I 

\ 


SCALAftS AMJ ARRAVo JL »J fcKPl-itIT DeFlWiTIUM - •**” MUT HcF£.HcUtfcJ JL 


*LL» t =»Nt!,OiJ if <‘«ii 

*H 6 =4iK 7 »NUMHP 

•50007 13 «SU005 14 •50005 

• S0003 ZO .^OOOif 21 .50001 

• 10006 25 •10005 26 .10004 

■10002 32 •lOuUl 35 *1 

♦120 37 40 


3 

fUl 

4 

♦ Nul-tlCI4 

6 

10 

3 

11 

16 

12 

IS 

•500U4 

16 

♦il56 

17 

22 

•soooo 

25 

•1C007 

24 

z'i 

.10003 

30 

*u 

31 

34 

•10000 

35 

TNA5S 

36 


T£HPaRARlE5 
•A0016 1/f 

IHL i NO fcftROMS. JETtCTEO X 






LA3ritL-FCiH 


FtJWTRAN VgaHtOiiii /Kl i*d-MAJr-7d 


03740 

OaYSO 

U3760 

08770 

08780 

08790 

0&800 

08810 

csaao 

oaajo 

08640 
086S0 
03660 
0887C 
08880 
08690 
03900 
08910 
08920 
08930 
0394 0 
03980 
08960 
0B970 
08960 
08990 
09000 
09010 
09020 
09080 
09040 
09080 
09060 
09 070 
09080 
09090 
C9100 
09110 
09120 
09130 
09140 
C91&0 
J9160 
09170 
09180 
09190 
09200 
09210 
09220 
092JO 
0 9240 
09260 
C92o0 
09270 
09280 
09290 


SiJbftijJTiNc: INPOTJ { 10 1 X • Y t 2 » T ) 
lMPt.lClT RtriU+ai A-Hta*21 
, CAUI_EO 8Y MrtlN 

COHrt«N 26XTRA/ WaDcAiNldtlFlcU 14i 

COMMUW /ELPAR/ NPAR Cl 4 Ji »1'(JMNP»HbAN0» NcLTYP tNl .N2» N3 *N4 • N5 oMTCr i f<dU 
dimension XCNUMNPlfYCNJHNPltZCNUMNPltltUNOHNPsO) ■TCNUHNPJ *IPRC161 t 
X ti(91 iA(9J»i. 11(9>*{.12{91 »C13(9l iC21 (9}>C22£9iiC23t91aV.31[ 

X 91 «C32(9} >C33(9J *844( 9) «a45t9) .864C91 iBSSlQl 4EU 191 »£121 

A 91 >L13 19) I E2U91 «E22( 91 « E23 ( 9 ) » 831 19 1 • E321 91 i £33( 91 , u44t 9] kJAsI 
X "1 p 6641 9) 1 u 65 ( 9 1 (HS {9} • RA(8 ■ 81 • RB Ca as 1 «RC £S«S 1 j EE (3 *8 J • 

_ X CtSiSl 

E 

: special NQDt; card fla>js 

E IT = coordinate system type ICC 1« ANV NODE CARuV 

: EQp » CARTES IAN 

» EW.C, CYCINDRICAL idlTH uEnENATION in CARTESIAN 

Z EQ*Di CVCINDRICAC trifrt GENERATION IN CVI-INDRICAL 

E IPK = PRINT SUPPRESS lDN ri.AU CCC 6> CARD FOR NUDE 1 LNLYI 

: EQ. » normal- printing 

: EQ.A. SUPPfttSS second PRINT IN5 OF NuDAL ARRAY DATA 

: £U*8o SUPPRESS PHlNTlNu OF iD-ARRAV 

• ED*C, BOTH 4A4 AND 

^ DATA IPRCXlrl a INA tlNBal HC*1H0/ 

IPR = IPRCi I 1 
NOME=NPAR( 21 
NUHL.SMPAR131 
NUHFX=NPAR(4 ) 

READ THlL.rt.NESS AND fcuASTIL PftCiPEKTl&S Dr THE LAVEkS, THEN DO THE 
Z TRANSFORMATION 

WRlTClOiS) ' . 

00 100 L=1.NUML 

read! 6«6l Aid *H( Lf *(«1 1 t L 1 • C 1 .cl 1.1 aCl S ( U) aCx^ZlL 1 vC231Ll »C33 (L 1 • 

X B44{L 1 fHASlCl >85S(C1 , 

C2l 14_)=C121LJ 
CjI CL l=ClJ(u) 

C42(C)~C23li-l 

ttS4{d=843(i.) 

C 

OU 96 1-1*6 
DO 9u J=l*6 
RAtI«Jl-O.DJ 

HOC I *J}=0.Da • 

RCl I . Jj=:O.DO 
L. 1 1 * J 1 = J *00 
95 Eeti ,31^10. DO 
C(1 tl)=LlUc} 

(.11 *21=C12(L) 

CM«Si=ClufUl 
CC2. D-1.21 {LI 



INPUTj 1_ASHi"u.FQ« HURTf^AN V.aAtoiiX) /Ki iia->JAr-7t» ljd:3tt PAG£; i-l 

I 093U0 > 

093iC Ct*i« JJ-Ciljtt-J 

oyz^a n=C3i<t-J 

09330 tl3*aJ=C3i(l_l ' 

09340 C13»3J=C33fu) 

093SQ Ct4^*4i=U44I£-J 

0936C C( 4»aj-rf4ali_J 

09370 CIS* 4)=Qb4tl-) 

09330 Cl5,b)=05SIt.) 

09390 KAU ■ l)=0CaS£Aft,) 

0940a RAIl .^J^USINIAICJ J4DSlflI AlCU ' 

09410 RAii i3j=“tJi>iNCa.4AtL} J 

0943a HAia,n=RAU»2) 

094^0 RAta,3)=RAU J 

09440 RAia,3J=“RAl 1 ,3* 

09450 RAI3.1 J=RAia»3a/-t:. 

, ' 09*ei0 HA(3pa)^RACi »3J/3t 




ItiPUTJ UiStiEj-mFiiH 


f-CiHlHAH /Kl 


iZiHa t*AusL £ 




HAlA*^)=aLQSiA£l^n 

«Af£i,4J=»ftA<4»53 

oa «6 i=if>s ; 

DU J=1 *iJ 
DO ^0 }4?=l»S 

RCliiiJ}^KCLl*Ji*ttAllvA7*CiK»J} 

HU tt 9X i~ftAti.n l 

ft0iz.ni-itAit*zi 

ni 

RUiSsli~HAil,M 

flS{ 

fibt 3 *JJ=«A( 3 * 5 J ■ 

3B£At4J=^JA{4»4) 

HtS{&t4i=riAiArUf 
HUi59£i~HAi4m 
CO 57 l=l9b 
DO 97 J^l*S 
DD V7 A~lwU 

E£iX*Jj=cfc£X*JJ+Rti£*A)*fllti£K*JJ 

Clt£4.)?=£Etl*U 

feXSti-J-E££lj>£) 

£13£i->=E££l»3X 

t2i££)=EEXiJ*XJ 

£.23tt.}=;£E£2» 37 
fc31tl.J-£K£3#l) 

£32 £t. 7 '£££ 3 * 27 
£ 33 («. 7 =££I 3 > 3 ) 

644£L7^Et£4*41 
0*fbtL.i=e.t.i4*Sl 
GS4££.7-££££«47 
CJibt 1 7 =£££&» i 7 

»ftIT£t6»771-#ri£l-J *£11££J »£lc£t.3 >££3fi.l fc£££JL7»l::22£jL7«£23£U7 Xa.7 • 
1£32£jlJ j£33£L:7«u-i4£LI«u4d(L.7 »Obb[L,7 

HSil3=rtUI 

IF .£7<i/4t. #£Q. i7 GB TB »ii 
OU lOl L.=2*«Urfi. 

£i=L-l 

Rt.mUD tb 

kHlTut lb) I XHlLJ *c i I U. J *fcl 2 £l -7 7 »£^ 1 £U I *c 224 l 7 

X i _7 <£ 3 i ti-i »£32 ti - 7 #t 34 ££l #*i 4 A£l -7 **» 4 Si jt. 7 >GS 4 ££1 »wSJ>£ 1 - I } * 1 -— I • 

4 ijltMLJ 

tiUAil irit. tlOUAL. COUHUittASS^ *iH IbT t-AV*,-* A»sl Trtc#< ^OHChAXC tut. 
poims Wi n(£ Tu'F tlF X3T LAYcft 

i7U 2Bti 





J o 


[— ] 




Lj-^ t..^J CZ-Hl c 


czj nzn 


INPUTJ 


i-ASHfc>.-*t-DR 




I2i3a PAGu 2-1 


1". 

10021 


T( Nl-O.OO 


10 022 

200 

R&AJ lSn300J liAUlTYiU.yi.UJ 


10025 

300 

paPMATU3.3F7.4i 


10050 

C 


1 

10650 

c 

K2A3 NUDAa paiNTS OHGREh OF FREEJUM lO 4V AbSuMlMo THEY 


10730 

c 



1074C 


DU IC4 I=1«NUHNP 


10750 


DD 104 J=sl,6 

i 

10760 

_104 

iDi l . JJ = 0 


10761 

• - %» 



10770 

■■ c 

READ ANV F1X2D NODE MUMo2ft A.N.j DIRECT laM 


10771 

■C 


’ 

10850 


IF {.NUMFX ,£Q. 0) GO TD 103 


10860 


DO 105 i=ifNlJHFX 


10870 

105 

READ {5(1061 Nt UDiNTOid J~l«3i 


. loaao 

lOd 

WRITE 16(2000} 

! 

10890 


write (6»2003i 


10900 


WRITE t6»200ii 


10910 

"C ■ 

- ■ • ■ 


10920 


DO 107 N=l*NOMHP 

■ 

S093Q 

137 

WRir5{6»2002) M( llDlNi 11 .1-1 .6) «A{Nl .YiNi .2tN> 


10940 

C 



10950 

c- 

generate the NOMdER DF cOUATiaM AND REDcFlNE. ID 


10951 

c 

FOR A SINGLE LAYER HRD5LEM 

) . 

10952 

c 

FOR MOLTI-LAVEH the PROGRAM SHOULD pc CNANuEO 


10960 

c 



£1000 


NEa=0 


IlOQl 


NUMfJl=NUHNP-l 

) 

§1002 


DO 6C N»1.NUHM1»2 


11003 


Ni=NFt 


1 1004 


DO bi J=1 .2 


11005 


1PU0(N«J)^11 62(63.61 

'• 

11006 

a2 

NEQ=NE041 


1 1007 


1D1N«31=:NED 


11006 


GO TO 61 


11009 

Da 

I0(N»3)=D 

> 

noic 

61 

CONTINUE 


1 1020 


DO 64 J~l »2 


11021 


IF UO( Nl .Jl-1 1 6bt66»b4 


11022 

o5 

NEQ=iNEm-l 

J 

11023 


1D(N1«J1=NE0 


11024 


uU TO 64 


1 1025 

ao 

1DIN1.J}=U 


11026 

04 

CUNTlNUt 

1 

11030 


lFUDiNT3)-tJ 67.66(60 


11031 

07 

NtQ=NE0+l , 


11032 


IO(Nt3}-NcQ 


1 1033 


lDiNlT3l=NEU 

) 

11034 


GO TO 60 


1103b 

DO 

lDtN(3l=0 , 


11036 


ID(Nl»3J=0 


a 1046 

Ui) 

CONTINUE 

1 

11470 

V 



114S0 

2 

PRINI MASTER IN06A ARRAY 


1 1490 

i. 



11500 


IF {IPR «EU> IPRClUi .URslPR .EO. iPHci411 3li TO 75 


ARli l~kCB 


rtS 


I' 




pi. 




l.-_L 




A J rVffri'f t 


' *’ ■>'-•* ■! - iatwififihnii jv W rii Ti i‘)b "rt il liMrli’i- 





r~] . 


INPUTJ 


11510 
1 1520 
1 1530 
11540 
1I5SQ 
11560 
11 570 
iisao 

11590 
11600 
11610 
11620 
1 1630 
11640 
11650 
11660 
11670 
11650 
11690 
11700 
11710 
11720 
11730 
11740 
11750 
11760 
11770 
11 7B0 
1 1790 
1 IBOO 
11810 
11620 
1 1830 
X1B40 
11850 
11860 


Ft^HTRAM V.541o^ll /»tl 23-HAV-7U 


WRITE Cb »2004) (N» UU(K. U 1 1=:1 «6 1 • t<=l « NUMNP 1 
76 continue 

1F1HL0£X«£Q.01 OJ TO 74 
data porthole SAVE 

WRITE INTaj e tliJCN. Ij pl=l»6l *N=1 ,NUHNP» 
write 1NT8) (K(M»N=1 sNORNPl 
UR1TI-CNT61 { Y(N1 •NsltNUKNPi 
WRITE tNTaj (2(N1 •»=! yNUHNPl 
KRITc INTO} CT(N) «N=l*NOHNP) 

ENDFILE NT« 

REHiINU H 
WRITE LZJ 111 

[ RETURN 

74 CUNT INUE 
REWINO a 
WRITE (81 JO 

* RETURN 

5 F0RHAT(y'y'19rtlELAariC PROPERT1E5J 

6 F0RHAT{F5altF3. 1.9F6»0J 

7 FaRHAr(lj>.iX*F4.1>2X>13F9.aj 
105 F0RKAr(lj«31 11 

lOsJU format (A1 >16 iAl >P10»1 >FlU«3>FiO«l *16>Fa0.i ) 

2000 FORMAT 1//23H1NUDAE POINT INPUT OATA J 

2001 FORMAT {5HONOOE 3X 24riblUUNUARV CUNOJTIJN COUbS llA 

X 2JHNU0AL POINT CuuROl.'<ATfcS / 7H NUMBER 2X INX 4X IhV 4X JrtZ Jx 
X 2HXX jX 2rlTY jX 2H2212X IHX 12X IriY i<iX Irii 12X IHT J 

2002 FDRMAUI5>615,3F13*61 

2003 FORMAT t//2lH wENpRaTEJ NOUAl JATA) 

2004 FORMAT ly/lTHlEOUATlUN NUMBERS/ 

X 35H N A Y 2 XX TY ZL /( 7161 J 
ENO 


12:38 PA6E 2-2 


COMMON ULOt^Ka 

/EXTHA/(+20> 
HUObX 4-0 


4-0 

NTd 

+ 1 

IFILl 

+2 

C+30) 

+0 

NUHNP 

+ 16 

MuANU 

+ 17 

+22 
+2 7 

N3 

+23 

N+ 

+24 


NElTY^ <-20 
Ns 1-25 


N1 +21 

MTUT +2o 


5U0PRU6RAM6 CALLEU 
DCOS. USlN* 





Lm CIZ3 czD ir3 [z:^ c-a czj cza 


tz=d ' tnZd L — J. LZZ 


INPUTJ LASrtEL^FOB 


FDRTRAM V.tiAto^iJ /K1 il4-HAV-/a 




PAuiC ii“d 


SCALAKS ANl> ArcRAYi i. WU dKHi-lCn DEFINITIUN - “ »** NJT KEFtHtiHC-tJ J, 


• S0C<!0 

1 

L23 2 

•I0Q12 

24 . 

• 1001 1 

25 

•lUOiO 

26 

E2E 

27 

2 SI 

4NliMMl 

52 

• 354 

63 

r 

7S 

■S0034 

76 

4N 77 

E33 

100 

•50033 

122 

rtS 

12^ 

■S0032 

14S 

-*A 146 

.50031 

14 7 

CX3 

160 

•Su030 

17*1 

H 

173 

cl2 21S 

♦ iPa 

237 

344 

240 

C31 

262 

EHd 

304 

4NUHt-l 326 

665 

J27 

034 

351 

y 

375 

czi 

374 

lO 416 

E13 

417 

*J 

441 

•50007 

442 

• soQoe> 

443 

646 444 

t.32 

466 

•3000b 

5l 0 

^44 

5U 

•S00D4 

S33 

RC 534 

.50003 

616 

E31 

617 

A 

64 1 

. soooa 

663 

*NortL 6o4 

•50001 

665 

• 5000 0 

666 

cu 

6o7 

♦ NOME 

rii 

EE 712 

635 

774 

•50017 

10 lo 

C22 

1017 

X 

1041 

• SOU lb 10 42 

•S0015 

1043 

fid 

1044 

•SQ014 

tl2o 

.10007 

1127 

621 1130 

.50013 

1152 

• 1003 6 

llSa 

•S0012 

1154 

• IC005 

1155 

• SOOil 1156 

•100C4 

1157 

4L.1 

1160 

•SOOlO 

1161 

.10003 

1162 

1163 . 

.10002 

1164 

C33 

1160 

•looai 

1207 

Q4S 

1210 

41 1232 

• IQOOO 

1233 

E32 

1234 

♦numfx 

1255 

C 

1257 

L12 1341 

•SQ027 

1363 

•50026 

1364 

RA 

1365 

■ SDOEb' 

1447 

Eli 1450 

. 30024 

1472 

IPHC 

l‘*73 

.30023 

ISOU 

•SC 022 

1501 

•S0C21 1502 







temporaries 















s s 

.A0016 

1626 

•aOOOO 1627 

.aoooi 

1630 





IMPUTJ i NO tBRORS OeTECTEl) ± 



V'*' 










MAIN* l-ASNtu FDRTRAN V.aAtlJ^U XKi l^:4a Pa(jc. I 



11B70 iJ&ftaarii^C J-AVSriL 

aiaeo iMpj_xcir HcALtettA-HrO--<!:> . 

11B90 C 

11900 C CAU-3 avbhL.str^c 

11910 C CALI.EO XiY tLTYPe 

11920 C 

11930 COMHON /(iAft/.A(99y9J 

119A0 COMaqwV'HLPAH/ NPAfUiA),NUt1NP,M£lAND»NEt.rr?flHltN3.N3.NA.N5fHnJT.N£U 

Jlgso CCHHGM /JUNp;/ lPADtSluI<>0}«N6iN7afVb«fil9tXdlO 

119fi0 CCHMON /EATRA/' NCiOEXtNTd.HlOSV.NTlO 

11970 C 

11980 Nl-1 

11990 IFINPAKU 1 >tU*CJ oG TO ^00 

12000 N2=Wl+6*HLrHMP 

12010 N34n2+NJANP 

12020 N4=N31-Nt;MNP 

12030 N5=N4+NUHNP 

12040 N6=F{Sl-RLrriNP 

IHOSO N7 =N6^NPAP13J 

12060 NB sN/X-NPARtBl 

12070 N9 -Na4NPAR(3J 

12080 M £j=N9+WAR{ 3} 

12090 HM=N10tNPAH(3> 

12100 CAUL, 2RHQR(HH} 

12110 C 

12120 call AYSHLlAlNl J >ACN21«A£M3J *At NAluNJHFXPJ 

12130 C 

12140 RETURN 

121S0 C 

12160 500 VtRlTE 16>20021 

12170 NuHc-NPAR(2) 

12180 i>6 800 HH=1>NUME 

12190 call STR5C lA(Nl) >A(NJ) ,Ncu,01 

12200 liRITE 16.20011 

12210 OU 600 U=!Lr*Lrt 

12220 CALL 3TRSC (A(NlJiA[N3I.NEQ.iJ 

12230 WRITE {6.3002J NHtU.Blw(lJ .5 ICl 2 1 tBlu { 3 } . S lul •. J . S 16{ 6j *S1 1>( ul . 

J2240 K 516C7J.S1S1 81 .Slt.191 

12250 C44» STRE55 PURTHULE 

<2260 IFlNlOBV.Ea.ll 

12270 AWRITE (NTIO) Nri .L .BIG 1 1 ) «S 13C .Blul 3 ) .5 1 u l4 1 « 6 1 6 £±>1 . 5 l 1 . 

12280 X 5Iut7) .516(61 iSIu{9) 

12290 800 CuNTINUt 

12300 fit T urn 

12310 L 

12320 20J1 FBRMAT t/J 

12330 2002 FURHAT (• IS HE CL ELEflENr BTRt^SC S«,y'7, 

12340 X * LlUMENT LQAU FXXl FXX2 FTYl FW2*. 

12350 A • FXVl FXY2 FYZi FY22 FX2*J 

12380 3002 FURMAT {2Id.9F12«51 

12370 END 




CGHHCN 8LUCA& 
76AH2 (+4/0 JOJ 


i » 


LAY5KL l_A5Htu'*f-0Fl 


A +13 

/EI-PAft/(+30J 
HPAB +0 

N2 +22 

NEC +27 

/JUNK/l+61 ) 
LT +0 

ri6 + 5 * 


/EXTHA/t+AJ 
HODEX +0 


cm 




dj 


PURTRAK V.itAl«2U /Kl 2J“«AY-7ii 


12 S ja 


pAuC I— I 


NJHNP 

+16 

MbANO 

+17 

WEuTYP 

+20 

N1 

+2i 


N3 

+23 

NA 


NS 

+25 

HTOT 

+26 


Lrt 

+l 

L 

+2 

IPAD 

+J 

51G 

+A ' 


NT 

+SS- 

NS 

+S6 

N& 

+57 

NIO . 

+60 

- 

MS 

+1 

NIOSV 

+2 

rtTIi) 

+2 





SU0PHPGRAH6 C,Al-L.nCt 

AV5HL ERROR STRSC * 

SCALARS AND ARRAYS i. **♦« NO EXPLICIT yEPINlTiUN - *'X“ NOT HEFLREnCEO i, 

+H« 1 *50001 2 iSOOCO 3 +NUHE + 

TEHPCRARIES 

■A001& SS *00000 Sb .UDOOl 57 .QQ002 60 .aOOOS 61 

CAYSHC X NO ERHOKS DETECTED ± 


c:d 




U- 



^ E3 ^^3 ^3 


AVSHL 


1^380 

ia390 

12400 

1241C 

12420 

12430 

12440 

12450 

12460 

12461 
12470 
12480 
124B1 
12490 
12500 
12S10 
12520 
12530 
12540 
12550 
1256C 
12570 
12SB0 
12590 
12600 
12610 
12611 
12620 
12621 
12622 
12623 
12630 
12640 
12650 
126b0 

12670 

12671 
12680 
126B1 

1269 0 

1 270 0 
12710 
12720 
E2730 
12740 
12750 
12760 
12770 
12780 
12790 
12BOO 
12810 
12820 
12B30 
12840 
12U50 


LASHcl aFUK 


ruHTKrtN V.5A(62l) /Kl 2a-Hrtr~78 


12:38 pAilc L 


SUbRUUTXiVt AVSHi-tlO«Xt V«Z »NUMNP1 
IHPUIClT ntAU#d{A-HjO-Z) 

1_0GICAU £uPkT»ELPCH«6ENPRT.6eMPtri»^icrPGrt 


CAUl_S CAL8A^ 
CALLEO BY 1.AYSMU 


DIMENSION XtNUMNP] tYlMUHrjP} tZlNUHNPJ ■ 1J1N0HMP>6) g 
X QX(999) •QV1999} tO£1999J >rZ(999> gElUO) gEl2tO) # 

X £22t9) *E2^(9) 1 53 1( 9 1 1 c32i ^ 1 g E36< 91 g044t9} gb45(9l f064l 9) , 

X 813(9) >E2H9) g055l 9) gTllOO) gilt 9) tn519) 

CUMHUN /UHluOD/t-*FX(999JgFFV(9&JJ gFrZ(99Jl 

cOUHQN /ECPAR/ NPARil4) gNNNNilgf4pANDgNc.-TYP9Nl ghlZi N3g«4* NS.HTOTi NfcU 
COHHON 7EH7L8i24J,M0,NS,St24,24)gPU4*4) .XMi 24) » ST ( l2 fTTtl2g4) 
CGHHON /JONK/ErtOLt4.41 « ll g Jl gXl *I2t J2 gX2»L «M >N g i j J gK« X2i Xl g Y2g Y1 . 

X DXtDY.DXYgOVXtDX2»HCU*Ri «R2gR3>a4tR5gH6*ft7»R6*H9gRlO , 

X MI M g HAX* NDl F f MT VPg MTVp£# R , u 1 i gM 1 g T 1 g T2, T3 , T4 , T5 1 T 6 

UDHMDH /tXTRA/ HDDtX,MT8 

COHMQM /CTl./ GR AVgHUPRr gEl_PCH.6eNPRTgGii:NPCH*6EFPt.H. 

I i WICugMlMBNOt 1P1.T 


N2=N1+6*NUHNP 

N3=N24NUMNP 

M4=N3+NUrtNP 

N5-N4+NUHMP 


CONTROL INFORMATION AND MEMBER DATA 


HTYPE-NPARUl 

IF (MTVPE *EQ. 14) Oil TO 2001 
read (5.2002) R 
f2 FORMAT tF13„3) 

U NUHE=:NPAR(2) 

NUHL=NPAR(3) 

MRITL t6t20aOl NUMEgNUMl-«R 


RbAJ THE THlCKNtSS AND ElABTIL PROPkRTISS FKUN TAPE 16 


REWIND 15 

READU5) ( UltL) gHS(L) gtl 1 (L) tbl2(t-l «E1 J ( i. J g c2 1 (u J g£22(L) .E23(L) g 
A E31(L) gE32(L) gEJ3lL) * G44tL 1 to46(L.J g G54(L) g655(LU fL=( g 

S NUMLl 

DATA PORTHOLE SAVE 
lF{MJDEX>EO<gl J 

XWRITE(NT8} Un(ul gHS(l-l g Li 1 { L 1 g£l2 (Ll t dl ^ t L J s L2 1 (L J «£22 (L) g u26 
A (L) gEBllL) tE32(L) gG44 IL 1 g G45(L ) ,Gb4(L) gLSOiL) 

K ItL-lgNUHL) 


EL(.MuNl uOAU MULriPLlEK^ 


NLAD (5gldU3) EMUL 
WRITE (6g20031 cMUL 

c**4 data pubthule save 
If (MODEA.EU* I ) 

X WRITE INTd) EMUL 





AYSHL. LASHui-.FUi^ FUHTkAM /KJ 2J-MAY-7b liiiJie PAIjE l-l 

12860 C 

12871 *KITtl lti.200t»l 

12872 C , 

1EB73 C euEMENT AmD uOAI> INFuRHAT ION 

12074 C 

12875 DC oOO 1=1«KUHNP 

iaa7U FFXtlJ=Q,00 

12880 FFY ( l)=:0«t30 

12881 bOO FF2li;=0*D0 

12890 100 RbAil {3>100bl M p 1 1 f J1 *K1 iL *HCD« HL 

12891 IF (HU sEQp 01 uO Tli lid 

12892 RhAJ (btlOObl QA ( 1 1 ) p <3Y C 11 1 taXlJl I tOY ( J 1 } i OX (K.1 } pUY IKU » 

12893 A aXll2}»UYU21.0XlJ2)paYLJ21«UAlA2)pQyiK2l, 

12894 A UAUUtFZ(Ill»QZ(31>«FA[J21 IU41A1) tFZlKll 

12960 115 CiO TO UOltlOZf 1U3.1041* MOD 



13110 

101 

A2=X ( J i J 



13120 


xi^xiin 



13130 


Y2=YtAl) 



13140 


Y1=Y1J1J 



13150 


JX-1 



13160 


jY=il 



13170 


3U= 1 



13180 


6b TO 105 



13330 

102 

X2~Xllli 


t 

13340 


Xl=XiJU 



13350 


Y2=yiKl> 

o o 


13360 


Yl=YUn 



13370 


JX=1 

s 

1 

13380 


JYs-i 

hd Q 


13390 


JC=-1 

o e 


13400 


GD TO lOa 

o ? 


13550 

103 

X2sX(il) 


1 

13S60 


Xl = X(Jl ) 



13570 


V2-VUI) 

^ 1 . 


13580 


Y1=Y(KU 



13S90 


jX=-1 

§ > 

> 

13600 


JY=-1 

fe Q 


13610 


JU=-1 

g td 


13620 


GO TC 10 3 


13770 

104 

Xl=XUl) 


i 

13780 


X^^XtJll 

£5i CO 


13790 


Y2=Y(Jil 



13600 


Y 1-1 n 



13810 


JA=" 1 


> 

13820 


JV = l 



13030 


JU4- 1 



13840 

165 

12 = 1 11^1 



13850 


JZ^Jltl 


» 

13860 


A2=Al+l 



13861 


JX=X4-A1 



13862 


DV=Y2-Y1 



13863 


OXY=DX40Y 


1 

13064 


i>YX=OY/OX 



1386b 


0X2=0X442 



13066 


IF IHTYPd *feU. 14) oU TO 160 



13870 


A»U = UA8biaX40Y«< IR+Htl.) )/6*) 



t 


»fT‘ 




! 

I 'CZ3 rm 553 CZ3 CZ3 cn c~] -CZZi 



AVSHl. LASHkU«hOK FcrtTkAh v'*i»Atuaj} /K1 aj-iMAY-ru 



13CS71 
J3S72 
t3S75 
13B?£ 
13B77 
13B7B 
l3aao 
t3BBl 
13882 
13803 
13884 
138S5 
13888 
1388 7 
138BU 
13689 
13898 
13B91 
1JB92 

13893 

13894 

13895 

13896 

13903 

13904 

13905 

13906 

13907 
1390B 

13909 

13910 

13911 

13912 

13913 

13914 

13915 

13916 

13917 

13918 

13919 

13920 
13930 
13940 
1394 1 
13950 
lJ9b0 
13970 
13980 
13990 
14000 
14010 
14020 
14030 
14040 
14050 
14060 


AAa=OMp5C5Y« (R+HtUU/U* ) 

GU Tu 161 

160 ABl=3Aa3tuX*iJV/o. J 

AR2=$DA3S ( aVfiJitU )/8« 1 

161 AR3=0AB6U1X4H1u)/3« } 

FF2C ll>=FFZlll)*Q2tll J*AB1 
FFXl U J=FFX(UJ+0X(1H4AR2 
FFYI UI=FFYUl J+QV(1U4AR3 
FFal JZ )=FF2( JlK-QZl J1 J+ARl 
FFX( JI)=FFaI J 1 1+OAlJl l^ARa 
FFY{ J1 l=FFVt J1 IFQYI Ji 1+AR3 
FFZtxn=iFF2t»U)+Q2lXU+AftJ 
FFMrtl J'-FF2iKl l+QX{^l l*rtH2 
FFYlXi J=FFY1K1 l+OYtM |4AR3 
FFXllZl^FFAllZH'QXtlal+ARZ 
FFYt J21=FFY( £2 1 +QY (121 4AR3 
FFX( J2J=FFX( J21+QXU214AR2 
FFV ( J21-FFY 1 32 J +QY ( J2) 4AR3 
FFX( K2J=FFX(K2I+QXlX214AR2 
FF V( K2 JaFFY 1 X2 M-QY (S2 J 4 AK3 
FFZl l2)siFFZUl 1 
FFZ1321-FF2(32» 

FFZ1X*;1-=FFZIX2J 

IF (MTYPti .60. IJJ 5U TU 107 

ai=M(t-> 

R2=l« 

R3=1./M(L> 

H4=H(l.J 

R6SHU.1/3. 

R7=rilLI/6. 

K8=H(i-)/3. 

B9=iJY 

R10=OY 

Tl-rUt.) 

r2=HlL|72. 

T3=HlL)/3- 

T4«Tl-2.*T^+T3 

T6=r2-T3 

T6=T1~T2 

3P=0 

GU TO 119 

107 R»-lRtri&(ull/n(U> 

R9=OYfl‘lR+HSt*»l J 
JP«1 

IF lU *20* 11' TU lOu 
Ll=U-l 

R1=1K*-«SIU1 1 )*HU.H-U.64HU.1+H<u J 
R2 =H+Hp(U 1 J+HCuixa, 

H3:=lR+HStl-l ) )/rm.l+0,5 
R4-MtCJ/(R+«6iUi J I 

R6=(R4-»f5(i.l ) /•*. 

R?-(R+MS{i-l 1 )*UU.)/6.+HU-14 M(i. 1/12. 
KOs^l R+HSU.1 ) )*HU-)/p,+H(L)*Htl-J/l2. 
M10-=OY»tR+nr,{H » )' 

GU TU lOJ 

lOo Rl~R4H(t.l40*a«ft{L}»H(LI 










AYSHL LASrtc*.-F0«5 f-URTHAN V.bAit>4i » /M ld:ja PAuu l--j 


14070 K2=R+fili,)/2* 

140BO aj=pr1/nti.) +0*b 

14090 R'4=rtti_j/K 

14100 Rb=R«rt(l-l/^«+ri(t-)*Hlt.>/4» 

14110 R7=R»ntL.}yb«+mL.)4HtUl/12. 

14120 n3sR4^n[C-)/3.-fHlU)«Hlui/IZ. 

14130 

14160 109 Tl=Ji-Uull.04RH} 

14170 T2=l *0-1.0/R44P1.00L1.U-I-R41 

14270 ro^O.b-^l ,0/H4+IDLUG<1 *0+Tl4U/K4/Ft4 

14320 T*.=T S-*.3»T24T3 

14330 TS-I2-T3 

14340 Tj&^TI'*T<. 

14341 C 

14342 C 1. f=3R« 5TIFFNEBS ANU STRB 69 HATRICcS 

14343 C 

14350 119 5( li U"0YX4(G44U,)4R3*DX2/12.+£lUul*Ra/2* J 

14360 St i* 2J"0YX*t (G45tL)+C»54U.) )#R3*DX2/24.+tei3tL>+£3Ul-J J*KU74. 1 

14370 Si !• j3=DYX/2.4(G44tUl*R340X£/ l2.-e.ll ti-)*R63+-{ElJll.l4E31iuJ >/l2. 

14380 X 4HU.14JC 

14390 Hi I, 4)=JYX/4#*tlG45lL.H-G54ti_J }«R5*JJi_/'l2.“(cl3il.H-ci3lCl_) J4ReJ>t£ 

14400 A 12ti.}4E21tUn4HU.)7i£. VJC 

14410 51 1, 5)=JYX*G44tL(*R340X2/24.-i£iJ4t.) lUJ 1*M(U)X12.4JC 

14420 5t 1> ol-UYX4{04blLl-t-G54(l.l}4R54t>X2/43.-(£12{i-J+b21(l.J J4H(U/1^« 

14430 X ♦JC 

14440 51 1» 7)=3VX/2**l-G44tLJ4tt34DXa/6.+£lU_)»R7> 

1445C &l I. ajs0YX4t-tG4btL)+G54tt.)l/R440X2/24.+ (tl3iLH-&JltL.J J*B7/4.> 

14460 St 1. 9J=DVX/2^*(-G44tL)4R2*i>X2/12.-tiU(t-14R7J + t£13tl.)+e3Ul,ii* 

1447Q X nli*)/24«»JC 

14480 Si l.lO|sOVX*t-{G45Ji-J+G54iL))/Fl4*DX2/45*-lsl3tU+EJltU)/4**Rr) 

14490 X tCbl2iUl+E21 U,] )4HtLlr JC/ 24 . 

14S00 31 * • li>--DX7*G44ti,J4KJ/24.-teiotl-J*£5l4t,» )4rttLl/24,4JC 

14510 31 1 f :2}=-DXY4tG4S(i-)-fGb4tt.n/R4745.-(£l^[i,)4i:2lU.U4JC4htU)/24. 

14520 Si 1 •l31=0Y43Y/l2#*lG44tl.J4R2-JPYt£t2tuJ+£2Ht,J14rttl-J/4-l 

14530 Si 1 tl4i=(0X*JX4tG4&{l-}-l-G54Cull-0Y4JV»2«4544£U>4R2-(£121l.J4E21 (l_l ] 

14540 X *JP»HtLI4DY*JY/2, )/24* 

14550 si A •l3j^''JP4DY4J74t£12il-]i'e2£tLn«rit.^l/4a.-0X« (Ci4btL>4Gb4tL.n/24, 

14551 X *JX 

14560 bllf l£>>-5UilJl 

14S70 5( 1 * 17J=5( !■ 141 

14580 Stl*ls)=Slltl5) 

14590 5t t.1 2] =jVX*lG55ti.»4Tl + Rb*RS*aX2/l2* jtLJ*rt6/2.1 

14600 51 2. 3j-UYX4i tG43ll-J+Gb4lL:J )*rtb»OX2/<»e« -ti.l3tl-J + b3l ii-n4R3/4. J+eJ 

14610 X 3tL)4Htl.)/6»4JC 

< 14620 

14630 Si 2* •*)-OYX/2.4tG5S(l.l4Rb4RS*ri*UX2/12«-£JjlL.)4Rti>+iE23tl-l4H32tl-l 

14640 X »*Htl-)/12.4JC 

14650 51 2* s]=DYX4(G45lL.H-Gb4(l.l}4R:i4DX2/4S.-e^3tt,l*Hll.}/6..4JC 

* 14660 Si 2. o)=DXY4GSbtl.l4RS4flS4Tl/24.-t£2-*li.)+fij2ti-J 

14670 3l 2. /)«JVX4l-tG45lLl+G54M-))*R54oX2224.+ (i*lJ<L}*'E3l(I-)i4R7/4-J 

14680 Si 2* 3»=I5YX/2*+l-G5Sil_)*RS4Tl/K4*0X2/o* +E33tl-14H71 

14690 Si 2. ifJ=0YX74,41-t645iL} + Gb4tl_l 14rt»*2X2212.-t£131l_J+£3lil_J I4R7J + 

> 14700 X njJtU) 

147X0 S t 2*i0 l = 0YX/2-4{-G55iLl4RS4Tl/H4*Ux2/l2*-i:j3lL J*ii7j4't£23tl.f 

14720 ^ +£.a2(i-) J*HiI-17'24.»JC: 

14730 St I i l=-OXY4tG4blt.J4Gb4il_) 1+R5/48.^EJJ tU>»rtiUl/12,4JC 





A.ysHL iiA3rifct-*Fafi 


FURTBAN V.SA(62U /KI 


PAae 


Se 2.i2)--DAV+li&5tL|^Ria*ri/H4/2A.-t£<lJl_)+EoitLH+Ktl.j/i!4.*JC 
£4 (iS‘^StL>FGSAlJ_l l#Rb-(Ea3tLJ i-u32(ul J/2. J^rttuJ 

i *jp 

sc 2.14)=-^t>¥*JVTi C043CL.>+G5^tU) ) «R5/24* 4- f ( U )4C^2U-) >/4S, ) C U j 

t «JP4-0X«JX*CrS5Cl_)4T14R3/l^« 

St a^lSJ^-OX^JX^GbSlD+Ti^^KS/lE.-DY^JY^^tE^JH.* f E32tL) >*HtLjy4B. 
t ■ ^JP 

St2»|6i=GC2.13} 

St2ii7J=St2» 14J 
5t2.ia)=S{2^151 

St 3i 3)=JYX*tu44{LJ4B31«0X2/ 12 * +C i It i-)*R d/ 2. i •• l£l J t U) *-t jX I*-) l*rl(Ll 
t /G« *JC+t33(L. J4T4/OVX/2* 

St 3. 4j=3¥X/4*»t tG4SU.)4-S54ti.X J «Ra*3X2/6t tU J +E41 IL) X*Rbl-t£ 
t X2tt.]i-£21 ti-J+3«*£33 lU }«Mtc. )/X2«4JC4t£aatUXFE32{EH 

< ♦T4/OYX/4* 

S{ 3. 5X=DXV»a44(LJ*B3/24* + <Ei3(i-)+c31ti.3i4mL3/l2,*JC-E33(l.l*T*7 

<. OVK/a. 

St 3> 6)-3XY4f G4atLj>GS4tEX J4RS/4a*4-tbX2tE3+£2l CL) J«Htl-)7i2.4jC-lE 
t 23lL)+E3a(Ej X*r4/DYX/4= 

St 3i ^X-DYX/2.4tt-G44tl_l*R3*DX2/t2.-eiltL|*R7J+tE13tU*+E31tJ.i* 

=l 4M( t-X/a^.tJC 

St 3, e J =3Yx74.4t-t G4btk-I+ub4li. X )/R44UX2/l^*'-lc.l3ti-J + u3i(EH4R7X + 

4 Ej>3(t-)4HtLl7X2»«Jt. 

St 3* ^X^OYXbtRt X tU)«:R7--G44ti.)4R34L)X2/P* }/2. -( E13 (X.I t-EaX tJLJX^rituI 
4 /'12.»JC+E33tE)*Tb/DYX/2, 

St 3<10X=DYX4t-tG45tL.)4G54(LlX7K44L)X2/’24.4 tEtS IEJ4-E3X lL:n4fl7/4. X- 
4 l2»4E3J(L)+tl2lEi*E2l(L> »4Ht L) ♦ JC/44.+ ( E2dtE i +E32tL X X+TS/ 

4 3YX/4. 

St 3*1 X X=-DXY#i»44tLX*R3/24-+ tEl3tl.J«-£3UUJ J*HtEX 724.<' JL-£33 lU *IS 
4 /DYX/E. 

SI 3 ,l4:)--aXY*lG45tLJ+aS4luX i/R474d,+ CEi^t uX+cEUl.) )4Hll.)+3C./24. 

4 -tl;234El+£32tLi J+T&/LIYX/4. 

St 3.13X=o44(u)*R24£0Y*JY/X2.+JP*ttt2tLX+-i:2 ttuJ X *H lu J40V4JY/46.- < 

4 c.2^tEH-Ej2tE) )*JP4TS4.JK*JX/24« 

St 3 tl4X=3V4 JY4{-2*4G44lL)4R2+JP4tEX2ti.X X^EdU.. X l + Kt UJ >2. X 724 . *-OX 
4 4JX4tG«b(EJ4GS4tEJ-tE23(t.X4-E32tEJ<X4JP4TbX724« 

St 3* lEX-JP-^OY^JY^lElEluX+EEKi.) I *Ht i-X /4 3, +3X* JX* t-tG4i>tEJ +Ga4 ti_I 
4 X-tE23tLX*£32tLX J*JP*T6J/24, 

S(3i XbX^Std. Xo) 

St3#X7X=St3»X4J 
St3, 10>-Sl3sl5J 

St 4» 4X =aYX*l GaSt i_)*h3*R5*T 1*3X2/ X2.+£33tLX4bti/ii,)-tE23CX.l+£*»2ti. X 
4 X*rttL X/6.*JCi^e22tU) *T4/0YX/2* 

St 4 » 3X=UXY*to4S(Ej4Ga4tLJ J*H3/4a. + £33lLX*rtli-X/o-*JC-t£2JtEJ+£32t 

4 i.n*T4/3YX/4. 

St 4 » oX=3XY*oaat L)*R5+rtb*ri/24 • 4-tE23 (l. X 4E32 li. 1 X *Ht L X / X c.*4Jo-c22t>- 
4 >*T4/0YX/2 i, 

St 4» 7)-3Vx/4.*(-(ii4bCL14G54tL) J*rt5*3X2/X2.-iE131i_X+£34t£XX*R71 + t 
t dl2<Ll+E2Xtl.jX«H(El*JC/24« 

SI 4 . bX=3VA'M-GSSti-X»Rb*Tl/H4*aX2/l2.-333il.X*H7X/2.+ lt2otJ_X+c,J2ti_ 
4 X l*HtLX/24 »*3C 

St 4t V}=sJVX*t“ tS4btu)+tjb4lbX )*Ka*3X2/24.t- td la tLX +E41 tL-I X*K7/4. J-t 
4 d«*ti33 IL:X+E12li_) t-ti2Xlt.l 1 *H(i.}*JC/24.f tE23tLX-i-E3^1EJ Xtlb/ 

4 3YX/4. 

St *t#I0X-=,jYX/2.»l-Ubbti-)*>{5*IX/R4*UX2/b. +t33lL)*R7X-tE23t*- X + 11321LX 
4 J*(iidX/l 2 .*JLi-E 22 tL J*rS/DYX/2. 


St 3«XXX 

4 

St 3 >12) 
t 

S{ 3«13| 

4 

^St bfl'*) 

St 3«lbX 
4 

S(3i XbX= 
St3#X7X= 
St3,10>- 
St 4* 4X 

St 4 > b X 

t 

St 4 > oX 
4 

St 4» 7) 
t 

b I 4 a bX 
4 

St 4t V) 


St 4tlCX 


t 


Cn 






CX] 


•■'n 


r:!] ' ci 


c 


c 


1 — 


j 









AySHL 

uASHci-.FDK 

r 

1S27S 

isaao 

S( 

4*113 


IS2S0 

A 



153D0 

5{ 

4*123 

1 

1531 a 

. X ■ 



35320 

5( 

4. 133 


15330 

X 


15340 

5( 

4*141 

1 

15350 

X 



15350 

St 

4 *153 


IS370 

A 



15350 

5(4*163- 

- 

153^^0 

5 14 f 1/3= 


15400 

S( 

4tl5J= 


15410 

SI 

5* 53 


15420 

51 

5* 6)' 

1 

15430 

5{ 

5. 73 


15440 

St 

5 « d): 


15450 

5( 

5* 93 


15450 

X 



35470 

St 

5*10) 


154B0 

X 



15490 

SI 

6.11 3 


15500 

5l 

5,143 

1 

15510 

SI 

5,133^ 


15520 

SI 

5 *14 3 


15530 

A 



15540 

S( 

5*153 

1 

155SQ 

St 

5, 16)= 


15550 

S( 

5*173= 


15570 

S 1 

5,163» 


i&sao 

5( 

6 ■ 6 3 • 

» 

15590 

SI 

6. 7 3 


15600 

5t 

6 * d 3 


IS&IO 

SI 

6* 9) 


15620 

X 


J 

15630 

S( 

0*103 


15640 

X 


15650 

5( 

6*113 


15660 

S( 

6*123 

1 

15670 

S( 

6*133 


15600 

St 

0*14 3 


15600 

X 



1570 0 

SI 

6*153 

) 

15710 

S( 

6 » lo 3 — 


15720 

S( 

6*1/3= 


137J0 

51 

6* ld3= 


15740 

St 

7. 73 

} 

15750 

51 

7, d) 


15760 

51 

7. 93 


15770 




lS7tJ0 

51 

7*10) 

i 

15790 

X 



I5BC0 

51 

7*11 3 


15810 

SI 

7*14 3 


1582 0 

S( 

7*1 j3 


rDAClifiK V.SA(u2U /lU 




i2ZJti 


PA on l—» 


bJ2(L) J«rb/JVX/4* 


1.3 ♦ rs/i* .'x/2. 


Il_ 3 • JP/4 U. -t 22 i 1,3 ♦ J P+ J X* J A* rt»/l 2 • 

:i3V*jY/ii4.* ( JP4 1 E23 lt.J tbJ2<<-3 3/2 *- tG4b {u3 +Ga4 i i. 3 3 VR&3 »nU, ) 
+0X4 JX / 1 2 « * C G 5b ( JL J 4 T 1 * Rb-J P* «i22 1 U 3 ♦ Tb 3 


+(i55tl.}*T14K3 3 


143 


. . - V V“ 


urx/2. 


*■103411-1 3 ♦T3/DYX/4. 


41o 52(L3 l^JP'OTli 3 


1212(1-3 +£2l(t-J J*Mti-3*JC/24. 


Lj4Tb/i3YX/2. 


>c4n( 


>t£:n4(u3 *T5/UYX/2« 


( 045 ( 1-3 -t-554( L3 )4ri5 


14.*JL 


3 41 1-3 t-t.41 IL3 J4i1<i.3/l4**J«- 


,*JC 


L. 




nun az2 nm 




r~~i 


c~n 




c~} cm r — 1 


AYSMt. i.ftSri£;,,FCiW 


FuRT(<AN V.bA(o21J /KI 


£j-rtAy-7ti 


PA&t 


15BJ0 
15840 
1S850 
15651 
15860 
15B70 
156B0 
15690 
1590 0 
15910 
15980 
15930 
15940 

isgso 
1596C 
1S970 
1S960 
15990 
16000 
16010 
16080 
16030 
16040 
16050 
16060 
16070 
16060 
16090 
16100 
16110 
16180 
16130 
16140 
l6lbC 
16160 
16170 
IblOC 
16190 
16H00 
16210 
1 6280 
162J0 
16240 
16250 
16260 
16270 
16280 
£6290 
16300 
16310 
16320 
16330 
Ui340 
16350 
16J60 
16370 


SI ) -CJX4 JAv 

X tu«Stl.|4'6b4ll-l )724. 

SI ?« 15)=OX«JX4cl04b(L.>+uS4tun/24i-6Y*JY«l612iU 1^621 iul 14H1L )/48 • 

A 4JP 

SI 7,16J=St7, 131 
5t7*17»=5£7.141 
S(7> 161=5(7.151 

S( 6. 81=OYX4C655(l.l«T17R4yR4*JX2/tt<±.FE33(i.l4Kb/2. J 

SI Q« 9>=UlfA/4.4( {045(i-)i)-GS4 (LJ )7R44 JX2/ 12 •- 1613 11. 1 FE3l (L) )*R6)i + 

X 233(i.l4ll(Ll/6c4JC 

St a*10)=OYX/2.t<(G5b(L|4Tl/R47-14vDx2/l2« -'633 U. iX^Ro 1 + ( 523(1- »4-£32(l. i 
A )4mi_)/l2,»3C 

St Q »i 1 J=i>XV* 1 o45 (l-> (t.H ( LJ *n(6 i/6 » *Ji. 

SI b el2}=UXy4655(LJ*Ti/R4/R4/24.'-(E23(^J i-L3£.tt,M*H(Li/12«* JL 
S( 6 .l3)=-JP*UV*JY/24.+ UE2b(L>+£32lLn*HlLJ/2. + (ti45(L3+a&4(LJ 1/R4 
A 4H(Lli 

St 8*14 »=«>Y*JY/24.*t (64b(LU-6S+(6 JJ/a4-JpjrHc<ijlLj-»-E32tLi J/2. J 
X *HtU-0K*dXCG5btLJ4Tl/R4/12* 

5( 8.l5j=DX+JX+G55tLj+Tl/ft4/12.-JP*0yTjY41£231LJ + E32(Lj i«mLl/46. 
btd.l6>=S(3.13) 

S(6« l7i=S(B* 141 
5(6.181=3(6.15) 

St V. 9 J=0XY4G44(1J 4R3/ la.-t-UYX* b 1 1 (uJ*R6/2 •- (£13(L ) +£al (Li J4H( Ll/O 
X *4jC-i-£33(L)#T3/2./0VX 

S t 9.101=OXV4{645(L)i-G54{U) l/R4/24«M£LalLl+Lal(LJ )4R6«DVX/4*’-(Ll2 
X (LJ+e2l(Ll+2.4£;33(Ll34H(6)/12.4JC*-4L23(Ll*E32lLJi*r3/4./ 

X OYX 

S( 9.11 )=UXV*G44(L)4R3/24.F(E13(L)+E31(_1 l«K(i-l/12o4JC-E33(Li4T3 
X /2./DYX 

b( 9.i<rl=UXV4«G4b(|.)+G54(Ll )/R4/4a«^(EiI2(i.li-£aI(Li }4H(L)/l2.4:JC 
X -(£43(L>-S-e32(Ll J4T3/DVX/4* 

Si s».131=0Y4JY/ia.*t-G44{Uj+R3+JP*:ttl2lw Jl-caKuJ ) /4« 14H(L)-OX4uX« t 
X cEblLlj-ESatLl 14JP4T2/24* 

St 9.l4l=dP*DY4JY746.4MtL141£i2(L)+E21(t. l4>4.«G4tf(L)*R3)-DX4‘JX4 
X t645(LM-6b4£LH't£23tL)+E32lL) 1* JP4T2J/24* 

St yilb)=JP*OY4 JY*H(LJ*(El2(Uf£21tLl J/48,+OX4JX4(G4b(LJl^Gb4£Ll 
X -JP4U£3iLj+h32{Ll)4T2»/24. 

St9.16}=St9. 13) 

3(9, I 71=3<9. 141 
5(9.161=3(9,15) 

SI 10 • lO)=dYX/2.« (GSblD^l l/R4/4f#aX£/D. «-£3 j( L )4R6 )— (E23(U1 +E32(L1 ) 
X 4H(L)/t),*JC+E22tLl*IJ/i>YX/2. 

S( 10 • 1 I )=OXr4 (G4b(L) -1-354 (Li ) /R4/4b.+tr3o( L > 4^Ui_ )/6.*JC— tE2btLl +£j2 1 
X L) 1 473/DYX/4. 

S( 10. i2)=i3XY4G5SlL14Tl/H4/R4/24*-H£2Jli,l -htoElLl )4HtL)/12.4JL-E2^tt. 
X )/2,/OYX*T3 

3(10*13) =OV* jr/24** t-t G4b(L 1 +G5 4 ( L) J /R4*- ) +E321L) 1/2 • 14HiL ) 

X -0X4JX4E22(L)4T2/12. 

5( 10 ,14)=3¥+JY/24.*tt623tLJ+d32(i_J J/a.Fl u4blL 1 l'Gb4 tU ) J/R4 }4H(Ll”iJX 
X *JX4(E22(L)4T24-G5S(L)«T1/R4J/12. 

S( 10,lb)=UY+JY/48*l‘tt2JtLJ+E32lLj )4irt( LJ *-0X*JX*GS5tL)41 1/R4/12. 

X -DX#JX*H23(L)*T2/13. 

5(10. 161=5(10, 13) 

5t 10*171=5(10.14) 

3(10.161=5(10,15) 

5( 11 .1 l)=0XY4G44tLJc:Ro/i2.+E33(L)VT3/iJYX/2. 


si 


p.hr 

I 

cr 


1 

^ cr 


C 


AVSHU 


1638U 

16390 

16AC0 

16410 

16420 

16430 

16440 

16450 

16460 

16470 

16471 
164Q0 
16490 
16500 
16510 
16520 
16530 
1654 0 
16SSQ 
16560 
16S7Q 
16560 
16590 
16600 
16610 
16620 
16630 
16640 
15650 
16660 
16670 
16680 
16690 
16700 
16710 
16720 
16730 
16740 
16750 
16760 
16770 
1678U 
16790 
16800 
16810 
166ZQ 

16830 

16831 

16840 

16841 

16842 

16843 

16844 
16840 

16846 

16847 


l_A3Hcu^F0rt FQUTRAN V^*6A£62U VkI 23-rt47-78 12J3b Hrtu£: i-7 


si 11 • l2)=aXV»iij4 6tl-l'*-6b4 il-l J/R4/*:4* + t£2s i*T-r/i*YA/4« 

SI ll»l3)J=-0Y*J7*644tt,)4R2/l2.*JP*6X»JX«<c23tuJ+e32U.M4T2/24. 

St 11 *14>=DV*JVi!<644Ii-J*aH/l2.-0X*JX*i J4blL) ♦654lUl’“(£2jtt-)+cJiitt.J J 
X 4JP4T21/24* 

S 1 1 1 1 1 6 j =DX«: J X* t G45 ll. 1 -i-GSA t u 1 -l-J P« i iT2 J 1 L. ) -tS lu 1 1 4 Tfi 1 / .<; 4 • 

SI 11 .161=Stlljil31 
S{ll«in=&ill«143 
Sill «16}=SU1»1 dJ 

SI 12,121-DXY*G5o(J_14Tl/R4/«4/l2«*fc22tl_J*Tj/OVA/2. 

Si lZ»13)=-0r4JY#(G4SIL*+Gb4tun4>rttj_J/S4/24.+i»X4JX#e22lL14T2/li:. 

X +JP 

SI 12»14)=D7*JV>l'lG46lt-J+G54ti_J l#Ml«-i/K4/24, 

X -[>X*JX¥l6abtul4Tl/34-JP*52bU.14T2)yiZ- 

St 12.1 3) =*DX+JX4iC5btI-J+Tl/R4-t-JP4t22t 1.1*181/12 . 

5(12.161^3(12.131 
5( 12«17)=3(12«141 
SU2.181-5tl2.lSl 

S (13 , 13 1= JP4DXV*e224l. I *T 1/48 <.-4-3 YX4G44 ( U) *R 1/6 . 

31 13sl4]=JP48XV*522(i.l*ri/9oa’’8YK«<344(6l *R 1/8 .-i- IG aSCL. 1 -t-6s<r 11. 1 1 4ri(i-i 
X J/16.+JC 

St 14 rib; “JP*DXV*EiT2tLl*Tl/96,-tG4&(l-l+5S4U.J )«mL.)/lb.*4C 
5( 13. lol=SU3.131 
SI 13.171=8(13.141 
3113.181=3(13.151 

S( 14.14l=JP«DXY«622(Ul*ri/48.-f‘3YX4G44(t-l*Kl/a.'-(64Sii.)-4-654ii.ll 
X 4H41-)/8*4JC-e-6SS(l,)4Ti/a./DYx 

31 14 .lb}=JP40X¥4E22(l-l 4Ti/96*-H G^bilLl ^o3«( 6} 1*H(L )/lo«*46''«»b( J-i 
X 4T1/8./DYX 

5(14.161=5(13.141 
S114,i7>=5ll4*l4J 
5(14.161=3114,151 

SUa,i51=JP*DXY*E22lL]l*71/4a<,-0-G55lU4Tl/a./iJVX 

5115.163=3(23,151- 

3(la.l7t=b(i4»lbl 

3(15.181=3(15,151 

5il6, 161=5(23, 131 

31 lo.l 71=S(13.24{ 

SUt>.261 = Stl3.tbl 
5(17,171=3(14,141 
3i 17.18}=3(14.151 
SUu. 181=5115*251 
05 111 jC=1.18 
08 Ml Xi-=l,ja 
121 StXC. JCl=5tJC.K.Cl 
f-U=i- 

OU =4B 1=1.9 
OO 349 3=1.18 
349 ET(L.J]=0.08 

IF 1 HGU .bU. 41 60 TO 334 
IF (MCO .to. 3l Gu fU 331 
IF (R60 .CO. a) GU TO 332 
Gu To 33J 
331 3X=-'OX 

R9=-R9 - 

Rl0=-Rlt3 
80 TO 3=3 


CZD CZD 





AYSHL LASHEL. *HtlR 


PuftTKAN V ,5 a (call aj-i^^Y-To 


i 

i 

f 

i 


! 

! 


16S4U ^^2 

16849 

16850 334 
16S51 

16859 333 

16860 
16870 
16880 
16890 
16900 
1691Q 
16312 
16920 
16930 
16940 
16950 
16960 
16970 
16980 
16982 
16990 
1700C 
17010 
17020 
17030 
17040 

17050 

17051 
I706Q 
17070 
17080 
17090 
17100 
1711 0 

17120 

17121 
17130 
17140 
17150 
17160 
17170 
I 7180 
1T190 
17191 
17200 
17210 
17220 
17230 
17240 
17250 
1 7260 
1 726 1 
17270 

17290 340 

17291 

1 7292 


DX=^-0X 
<5D TJ 33J 
H9=-RS 

mc=”aio 

5T( 1>11--611{L1/0X 

STl i«21=<-£13tL)2DX 

ST£i f3)=eiUul/DX-Ei31LjyBl0 

ST(ly4J=-£l*:(LJ-7R10+E131L>/DX 

STCA»5l=t3l3ii-l/'Rl0 

STU o6)=£l2{L.l/Ht0 

STtl »131=:JP»Ei«ltLJ/R/4. 

STll t.l4l=JP457U , 13J/2. 
ST(2»rj=3T( 1 ill 
ST(2p8)-8TU*2) 

S T ( 2 . 9 ) ■= EH { t. )/ DX-t 1 3 1 L) /il9 

STlailOi=~El2t»-J^R9+E13(LJ/OX 

ST12.M3-E131LJ/R9 

SIl2Bl2)~El2ti.0/R9 

ST12»1 j)=JP4E12t;Ll/44y{ft+HSlL» 1 

STl2il4)-JP*STl2. 131/2* 

STI 3 .1 1=-E2l (LIZ-OX 

ST 1 3 1 2i=-'E2ot J.I ^DX 

ST43 ♦31=£2l {L1/0X-E23CL1/R10 

ST 1 3 .4) =s-£22 tt.J /R1 0+e 25 11. 1/0 X 

ST43»51=E23{U)/R10 

ST13»6I=E2214-1/R10 

ST13 •13J=JP4E22a.I/4-/R 

3ri3>14)-ST{3*13)y'2* 

STt4.7)»ST13.11 

ST(4«8»-ST(^i2> 

SrH«9)'=E21U-)/l>X^£23tLl/K9 

sr 1 4 * 101 =-£22l4u 1 /R9+E2311.J/OX 

ST £ 4 ill l^ERSH-J/^RS 

STt4.121=£22U.l/R9 

ST 1 4 i 1 jI= 3P«£22(1.) /4-/( R+HSlL.1 1 

ST £ 4. 14)=ST( 4,131/2* 

ST(5«1 ]=-E31 (UJ/UX 

ST £5 ,21“-e33 £U J/DX 

STt5,31=£aULl/OX“£33l»_J/Uitr 

ST(5,4J=-t32St,»/Rt0*£33lLl/ilX 

SU5»SJ=C331J-J/R10 

ST(6,6)=to21l_J/KlO 

STl5.13»=Jp*E5211.1/4./rt 

STi5 ,141 =SHS, 131/2. 

ST16,7)=STC6, 11 

ST(6,0}:=ST(6,21 

STl6»9J=ii31 lU)/l>X-£33 lLjyR9 

ST 1 P e 10 1 is-fc j 2 1 1. 1 /H9+t33 1 J_1 /O X 

ST l6*l U=£3ilU}/R9 

STl 0,1 21 =LJ2U,)/R9 

ST lb,13>=3P«£32(£l/4./lH+HSU.M 

6 T ( 6 « i 4 1 — S T 1 6 ■ 1 P 1 /2 a 

OU 348 MJA=laU 
ST l8JA*l5l=ST£MJAa l4) 

Du 347 MJC=i ,0 
DU 347 HJD-I 7,18 
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AYSHL i.ASHfcL*{*UR FurtTHAN V.oAto^ii) /fU -ij“MAy-7«i J-V 


* 17?,$I3 347 ST{MJC*HJJ)=ST(MJC*i4} 

*r2S4 UU 34b 

17295 346 SrtHJE»i61=5Tt‘i|Jc.l3) 

17300 S7(7.U=-^iS4tU72.7rttt.> 

• I73I0 sri7p2J=“65Bli_J*89/'3.7rt(l.JVRJC 

17320 STC7*3J“5TI7,l>y2- 

17330 ST{7.4J-3T(7*2J/2* 

17340 STC7.5)=ST£7»3) 

I 17350 an7f01=br£7*4J 

17360 ST(7 + 7j=-STt7»U 

17370 ST£7,dl=u5b£l,)72,7H£t-l 

173B0 STt7p91-=3T(7»7>7a, 

' 17390 STt7ilCj=ST< 7.BJ72P 

17400 STf7*Il )-ST£7»yJ 

17410 STl7slZ)=STl r,10) 

I 17420 ST£7.13j=-GS4tl-J72*/l>X 

[ ' 17430 5T£7»I4J=*-SS5t«.372./R10*-6ii44LJ/2*/3X 

17440 STtr.3 5l=G5a£l.J/2.7RlC 

: 174S0 ST£7pl6»=STt7#i3j ‘ 

* S7460 STC7.17)=STt7»141 

17470 ST(7.1UJ=3T£7#131 

1748U 5TlO#l)="G&4£i-J72.Xo£l,J 

17490 ST£e*2J='-GS6£i,j72,7rt£l,J 

I7S00 57£a«3)^5r£dall/2» 

i ' J7510 5T£a.4J=3r£d«2)/2. 

i 17520 ST£0«5)=ST£a>3| 

( 17530 ST£a.6)~3T£B,4( 

17540. ST£0»7J--ST£ a.l J 

< 17550 ST£aiaj=G55£t.Hiftl072*7'l£l.J7ft9 

17560 ST£8p9)-5T£af7J/2. 

17570 ST<d»101=5T{a9aj72. 

irsao 5Ti6#n>=ST£a.9J 

17S90 STtaal2}-=5T£a.|01 

t 17600 ST£a»13;=-G5ACl_)72.70X 

17610 5Tfa,14J=-G55£L)72.^RS+GS4Cl.l72.70X 

» 17620 ST£a»161=055IM72.7ftiJ 

I 17630 , ST £aaI61=Sr£a.l3J 



17640 

ST(a«17J=ST{a«141 


17650 

ST£a*ld)^ST£a,l3> 


17660 

Sriy.lJ=-G44£l-)72.7MU.J 

J 

17670 

ST£9a2)=-u46ll_l»R*//p-«/Mlt.J7«lD 


17600 

5T£9,Jl=5T£9tl J72, 


17690 

{9.4l±=oT{9«2)72« 


17700 

5T£y,sj=srf 9*j| 

J 

17710 

ST £9*6J=ST19»4) 


17720 

ST l9»ri =-5T£9.i 1 


17730 

5 T £ y * ») =1*450. 1 72.7H£j_i 


17740 

ST £9.9)=£Tt9p7)72. 

1 

17750 

sTty»io?=5Tiy.e}72. 


17760 

SUV.J 1 J=ST£9.VJ 


17770 

51 ryTi2J=sr£9»j07 


17700 

ST£9.1JJ=*-G44£LJ/2.70X 

1 

17790 

SI £9.I4J = -b45£l-l72.7Kli>+G4tO-l72 


17300 

Si £y *15J =(*4 u£JLJ72 p/R10 


17010 

ST CV*16J=ST£9*lJJ 


17820 

ST £y*17};^T £9* 141 
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AYSHL l-ASHc_ .FUrt FORTRAN Y*6A£oiiii /^r 2^**MAT-/3 1 E; 


l?a30 3Ti9ilo)-S7l9tl:>| 

17640 C 

17650 C 2. aET ORAtf 1 TV LOAD rtArit ICES EO UAL ZtBu 

17660 C 

j J7670 00 450 I^lalS 

17660 DQ a50 C=l*4 

37690 330 pl£»Ll=0>00 

17900 DQ 375 i.= l»ia 

J 17910 375 XH1U>=0*30 

17920 DO 47q t_=l*4 

17930 DO 476 J=1.9 

17940 376 TTt3»ul=ai.OO 

17950 C * 

17960 C 4> FCRM eGCATIUM HATRiX AND CQ.4BUTE dANO UIDTH 

179/0 C 

17980 DO 400 i^=l>2 

' 17990 

i 18000 LHiL-^2)^I0t41iE) 

16010 i-NU.44)slO(tU«LJ 

f laoao LM(L*-6J = I 0 tI 2 .O 

t 16030 LM{ui-6) = I0{J2.E} 

18040 400 LMlL+10)^IDiK2«i-} 

' laOSO LH(lJ>-IO(lt«3) , 

18060 EH1141~1D(41 >31 

■ i 16070 UH(l51 = IDlia »31 

f ‘ 18080 LH(16J=I0(12,31 

f 18090 l,KU7>=IiJl4a.3» 

iSlOO LH( iai=2D£K2t3) 

: » IBllO NDH=24 

18120 ND=ia 

! 18130 NSss9 

18140 IF l.NOT. Ei^PHTl OD TU 408 

I ' laiso HRiTE i6«aoa7J ndbNS 

I 18160 IF (£t_PCH) MRITE (ll>20Ca) N*K3aMS 

i 18170 00 406 1.-1.N8 

\ 18180 IF (ELPCHl 4R1TE U 1 •2D09l(STCL*|Ci aX-l iNQl 

J 18190 406 ItRItTE (6>201 0 J 1 ST IL.Kl *K= 1 cNDJ 

r 16200 IF (£3_PCHi bRlTE U 1> 2089) ( (TT (L aK.1 tK= I » 4} * 1.-1 *NS I 

18210 «R1TE (6«20lf MlTTCi.t«O.K»i*4)»l-=liNS) 

16220 408 CALU CAE8AN (f48AND*N01F»LH»XHii3»P>N0.N0M*l4S) 

, i 10230 NTl=l 

18240 IF {HiNbND .EU* 0 1 uQ TQ 409 

18250 NT1=4 

18260 EHI1=l.MI 134-LNl«i)^l.»Ul3) 

) 18270 l.rtJl=I.HI31+LM<4J *LHtl4) 

16280 LHK.l=UH{S)+l_HU>)+iLM(15) 

18290 LH12=EM17)+I.H18)*LHC16» 

16300 1.MJ^^M19)+LH{10}+EM(17) 

) 16310 LHK.2=l_M(l U+LM( 12H-LrtCiaj 

18320 IF ( a.Mll«EHJl«t.HXi4|.Hl241.MJ2*I,Nr;2i .NE- O) 

18330 \1«R1TE (jiJ 11 *J1 •Kl«£2tJ2>K2 

18340 409 LilNllNUt 

} 10450 IF I IPLT .oT . U> 1«KXTE (13*4 ) i 1 ,41 , iv.1 • £2 > Jci 

16460 lF(HaDEX*tm>01 TU 410 

1B47C UATA FURTHULE SAVE 

18480 MRITE (NTd) M, 1 1 *4 1 ,K1 « 1 2*42 aK2,l.l 1 










AYSHL LAShC-.FUft (-URTKAH V*5Aiu2l) /K1 12i3et P/iGC i-B 


18390 liLi T*i 42U 

I 1Q40D 410 i-ONriNUE 

I 18410 )ifRlT£.CNn fi.=l»Ni31«CC8Tu.iAl ,i.=l,f4&Jt<=liNO}« 

I IB42D X ( tTTU.«j<l «X=1 f4) 

t 18430 C 

I 18440 C 4. CHECK FUR HURE ELEHEHTS 

I 18450 C 

I 18460 48^0 CQHriHLlE 

I 18470 HRlTE Cbf2004l Mu 1 1 • J1 tKl •X2« J 89 K 2 «l.i 1 iND IF 

' 18480 U=Ll 1 

18490 IF IM *Ea, NUHLJ RETJRM 

18S20 3U rc 100 

[ 18530 C 

i 18540 1003 FORHATl4FlQ«i) 

I IB55C I00t> FOBHATt4l&»JtlJ 

18551 1006 FflfiHATJ 18S^4.2 ) 

18560 2000 format t'lSHELl, ELEMEN T 

18570 < * NOHBEK OF SHELL EUEHENT&=* • 15i /. 

t 18580 X » NUMBER OFOIFFp LAYERS** . 15./ • 

> 10590 X • RAOlUS**.F13*J) 

^ 18600 2003 FORMAT L///25H ELEMENT lUAU MULTIPLIERS / .cOx *1HA« I4X« IHO. 14X. IHC* 

y- 18610 X L4X«iH3./6n X-OlR.4F10« 1/ faH Y-UiR.AFlO . 1 / 6U 2-UlR,4F10a 1/ 

[ 18620 X 6H T£MP»4F10.il . 

1B630 2004 FORMATtai5*I7> 

16640 2003 FaRMATi///47Hl « 11 J1 K1 32 K2 LAYER uAUO 1 

18650 2007 FORMAT I'ONuHaER OF DISPLACEMENTS PER slEHENT tNDJ =«*13./« 

! 18660 X *0NUHQ£R OF STRESSES PER ELEMENT twSi 

i 18670 X 'OELEMENT STRESS-DISPL MATRIX »J 

f 18680 2008 FL^tMAT A* ELEMENT • *12 »* ND=**i3*e N3=* .151 

18690 .jOOy FORMAT I C1P9E10.311 

18700 2010 format iX.tlH .1P10E13.41J 

18710 2011 FORMAT 1*0£UEMENT FIXED-NODE STRESSES *,2.Un .1P4L13.411 

IQ 720 END 
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XI 
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Yl 

+62 

OX 
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+0 6 

DXY 

+70 

DYX 

+ 72 

0X2 

+ 74 

MCO 

+76 

R1 

+ /7 

R2 

+101 

R3 
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H+ 

+ 106 

R5 

+ 107 
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+111 

R7 

+113 

«a 

+ 1 l3 


+ 117 
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+121 
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+I 20 
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NDIF 

+125 

HTTP 

+ 126 

MTYPd 

+lii7 

H 

+130 

Ltl 

+132 


HI 

+133 

T1 

+ 134 

ra 

+ 136 

73 

+140 

T4 

+142 


TS 

+1A4 

T6 

+146 
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✓EATRA/l+2) 










HODEX 

+0 

Nib 

+1 








/CTL/t+121 









t 

GRAY 

+0 

ElPHT 

+t 

fc-LPCH 

+J 

GtfiPR t 

+ >♦ 

ocNPlH 

+5 


peFPCH 


' TCG 

+7 

HlNbHu 

+10 

IPlT 

+1 1 
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CAL BAH 

DAOS. 
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) 

SCALARS 

and arrays 1 HO Explicit defimitioh 

- **A'* HOT KcPtREHCEO X 



•sDoao 

1 ■ 

»Jp 

2 

*AR3 

3 

OY 

6 

•10010 

3723 


E2Z 

3724 



4HJO 

3746 

4LHA2 

3747 

J4T 


J 

»SOOJ4 

3730 

E33 

3/51 

.50033 

3773 

h5 

3774 

•5QC32 

4016 


.50031 

4017 

#MJHWP 

40 20 

.50030 

40 21 

H 

4022 

OX 

4044 


E12 

7762 

4MJA 

10004 

4LHX1 

10005 

4HDM 

10006 

+JC 

10007 


Ek3 

lOOlO 

*MJE 

10032 

654 

1C033 

Y 

IC065 

+JY 

10056 

I 

*ARl 

10057 

10 

10061 

, E13 

10062 

4lHI2 

10 104 

.50007 

10109 


.SOOQG 

10106 

. 50005 

10107 

G44 

10113 

• 50304 

10162 

.50003 

30 133 


♦JX 

10134 

E31 

10133 

.50002 

10157 

4HJMl 

10160 

FZ 

10&61 


•50Q01 

14077 

•SOuOO 

14100 

♦NTl 

14101 

4NUME 

14102 

355 

14103 

i 

+LM11 

14125 

^Mu 

14126 , 

.50017 

14127 

X 

14130 

•SODlb 

14131 


AAH2 

14132 

.50015 

14134 

.50014 

14136 

.1UJ37 

14136 

E21 

14137 . 


.S0013 

l4lol 

. lOQOo 

14162 

.50012 

14163 

J w 

14164 

• JOOOb 

14165 


•seoii 

14166 

• 10004 

14167 

4LMJ2 

141 70 

• SOOl ij 

14171 

«L1 

14172 

1 

.10003 

14173 

. 10002 

14174 

•10001 

14176 

ut5 

14176 

QZ 

14220 


.10000 

20135 

L52 

20137 

.50027 

2C161 

•50326 

2Clo2 

.50025 

201d3 


Ell 

20164 

• 50024 

20 206 

+lH3I 

2020 7 

.60323 

20210 

«tcc 

20211 


.50022 

20212 

.50021 

20213 







1 

TEMPORARIES 
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20407 

> OOOOg 

2U410 

.00001 

20411 

•00002 

204 1 2 

•Q0003 
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cr3 Em crz3 1—3 CIS nzD c~] era ca 


err) err cm 
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f 

phintd 

LASHeL-FUtt FOUTrtAN y<bAlo2il /M 23-MAY— 78 12:38 


la/do 


SUBKUUTlNll PrtlNTO >r4t.ai>i 4JMNP«[.l.<Nu!.UL.K.«Nt=Q 

j 

18740 


1 HU) 


18750 


IMPLICIT Rfc:AL4atA-H*U-il 

r ■- 

18760 


real FR 

1 

16770 


logical £LMRTiELPCH*GEHPKT*u£NPCrt,J2Fpt-H 

1 

18780 

C 



18700 

c 

LALlEu 8V SULLQ« SOLCIuvRESPLi- 

[. 

16800 

c 


' • - 

laaio 

c 

HQsl static 5UL3T1UR ^ 

1 ■ 

18880 

c 

MC1=2 OYNAMIC SOluTIUR 

i -, ■ 

1B830 

c 

HQ=o KE3PUN3E GPCCTROH ANAi.Y;>13 

) '•■ 

18840 

C ■ 


) 

18850 


OIHEHSION I0iNJHNp«6J •iH61*F;h5d*LL) .F2£tlll«' 1 


18660 


CUHHOM /CTL/ URAV*fcLPRT*ELPCH»5EHPRr*aEHPCriiO£FPCH 


18870 

tl 


[ 

18880 

c 

REAO THE 10 ARRAY AND THEN HEAO THE GOi-OTION INTO CORE 

^ t 

18890 

c 


1 ■ 

18900 


REWIND a 

s 

18910 


1F( UbEU.BIMEWIND 7 


16920 


READ Caj ID 


16930 


REWIND NT 


18940 


IF IKU .£0. 3) OEFPCH=.FALSt» 

t . 

16950 


IF (HO *EQ> 1] GO TO 110 

' . 

18960 


1MMQ.NE.3) READ INT) FREQ 

: ■ '■ ■ 

18970 


IFtHQ-tQp3j READ (7) FRtQ 


18980 


IFlHO.EO.j) READ (NT) 


18990 


TPI=sa«OOU40ATANil aOOO) 


19000 


DU 100 N=1.NF 

\ ' 

19010 


FHtOfH)=FRE(itN)/TPl 


19020 

100 

CONTINUE 


19030 

110 

LONTlNUh 


19040 


NL=NS8 

1 

19050 


NS=HE+l-NEOtl 


19060 


00 120 N=I*NULUCK 


19070 


if-(Ha.E0.31 HLAOtri 


19080 


READ (NTl (CF(JtKl*3=NSiNEl«iv.=l«LL) 

t 

19090 


NS=NS-NEOd 


19100 

120 

NE=N£-NE08 


19110 

c 



19120 

c 

PRINT THE unrestrained NUDES * D l6Pi_AUE4ENT6/RuTAT ILnS 

t 

19130 

c 



19140 


00 199 L=1*U1. 


19150 


GU TLi U3U«14U>150)»HU 


19160 

1 JO 

WRITE (6t2000)L 

1 

19170 


WRITE {6«20JU) 


1918D 


GO Tu 160 


19190 

140 

FK=SN6L(FREu(L} ) 


19200 


WRITE lbi2U10t L«FR 

t 

19210 


write (6*20jOJ 


19220 


GD ru 160 


19230 

15 0 

FR-SNuLtFREaiL) ) 


19240 


IF tL .LT* LLl WRITE (6*^0^0) 1_.FR 

} 

19250 


IF (L .LD. LL) write (6t204O> 


L926Q 


WRITL (6<20jU) 


19270 

160 

00 IdO 1'=1>1 '|UMnP 

. 

19280 


IFLu^U 




PR£NTU LASHEi-.FQrt 


FlIRTHaN V.^AibHU /M 


-7u 


P*Vuu 1-1 


X9290 
19J00 
193X0 
19320 
19330 
19340 
19350 
19360 
19370 
19380 
19390 
19400 
19410 
19420 
19430 
19440 
19450 
19460 
19470 
19480 
19490 
19S0 0 
19510 
19520 
19530 
19540 
19550 
19560 
19570 
19500 
19590 
19600 
19610 


OU i70 J=li5 
Dt J)=0.0D0 
N=IOl 1 1 J ) 

IF IN -£0. 01 iiO TO 170 
IFiN.i-T.O mAP8=-N 
IP1N«1_T<0 )N=icHNAPuiJ ) 

1 Fi.G-1 
0( Jl^FlN.U) 

CONTiNUc 

IF 1 IFUG »tu. 01 uU ru luo 

ViRlIiI i tiJ 

IF CJcFPCri) WRITE (ilt^OoOl l*u.^ 

LUNTINUE 

IF (OEFPCH* «R1TE (11.2070) 
cam IMUu 

FORMAT (*15 T A T i C A N A L Y 5 1 5* • UOAO ClA5fc*»l6j 

FORMAT ( 

K»0Dl3PLACEMENT3/HUTATiaNS UF UNRESTRAINED NODtS'.yy. 

X* NuUL • i 2( l<iX, * A-* d 12X» » Y-* .tiJX. ) ,/• 

X* NUMoER* «3(3X. 'TRANSEATiaN* 1 ■■*(oX«*aarATlDN> ) } 

F6RMATC1M D A I- A N A i_ Y 3 1 S* * RuDE NUMDER'tib/t 
X • FRc:QUENCY=* .1PG13.6.* H2*.//. 

X * E1JENVECTOH5 NORMAUlZEU TO A JNlT MASS MATRIX*) 

FuRMATl* IRtSPOMEt SPECTRoH ANA 1-7615* 

X * RESPONSE FUR MODE* .15. ✓.« FrtfcJUtN^-Y=* .IPGl ^.a.* H2» > 

FORMAT ("IR ESPONSE 5PtCTriJM aNAI-YSI 3*.//, 

X * SauARE HOOT UF Trib SUM OF THc S0UAR&:, OF THE*. 

X » NOOAU deflections* ) 

FORMAT (io. lx. 1P6014.5) 

FORMAT (2I3.2X, 1P6D12.5) 
format 1/1 

RETURN 
END . . 
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GRAV l-J 
UEFPCh J-b 


ElPRT *-2 


cLPCH +3 


tlENprlT +4 


uENPCH +5 
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♦ TPl 

1 

FR 

3 

♦LL 

4 

«NE03 

5 

*NS 

6 

i 

4NbLULK 

7 ■ 

+N 

lu 

*K. 

11 

♦njmmp 

1^ 

XNbO 



♦ Nh 

13 

Lv» 

14 

lo 

15 

*3 

lo 

•SC006 

17 


•SO 005 

20 

O 

21 

•30004 

35 

4NS3 

35 

•30U03 

37 


4NE 

40 

•30002 

41 

•SOOOi 

42 

.3000 3 

43 

•NAP6 

44 

i 

. 1000 7 

4.3 

. lOCC^ 

40 

•10005 

47 

• 13004 

30 

*MU 

51 


• 10003 

52 

♦L 

53 

■ IU002 

54 

*1300 1 

55 

• 1. 

5b 


• 10000 

57 

rNuU 

oO 

h 

ul 

•NT 

6c 
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PR INTO LA5KEU.F0R 


PUKTRAN V.i»At621) /Ni /:J-MAir-7o 


s I 


TEHPCRARIES 
,A0016 260 

PRINTO i NO EHRDrtS OETECTCO ± 


Q O 

^ OQ 


i 


O 


iizri . nri i~i • [~i r~i tnzi rz i r~ i CZD 


PAiic i-c. 



L.ASH£i.*FOa 




12: 3d P^tj£ I 


i9620 

19630 

l964rO 

19650 

L9660 

19670 

19650 

19690 

19700 

19710 

19720 

19730 

19740 

19750 

19760 

19770 

19780 

19790 

19800 

19Q10 

19820 

19830 

19840 

19850 

19860 

19Q70 

19880 

19890 

19900 

19910 

19920 

19930 

19940 

199S0 

1996C 

19970 

19980 

19990 

20000 

20010 

20020 

20030 

20040 

2O0S0 

20060 

20070 

20080 

20090 

20100 

20110 

20120 

20130 

20140 

20150 

20160 

20170 


3UdKuJTlN5 ahSljL (A»B.i>IA4AtNcJ ,MA*N</«t4aL3us.f^bQb*N3V»Hl tNSI IFi 
I NrtED.NL.NRJ 

1RPI.1&1T REAt-4=aiA-rtiO-Z» 

■ CAU-ED er sobta 

UlRENblJN AlNAVJ •dlNAVJ (MAXAIMI) 

HHii 1 

MA2=MA - 2 
lFtHA2.EQ*0| HA2sl 
lNC=NEQd - 1 
NWA^NEOd^HA 
N78=1HA-2|/NEQ0 +1 
NE8=NTd*NEQa 
NEBT=NEa + NtQd 
NktV=UEQb4NV 
Nl»vy-NEdT*RV 

N1=NL 

N2=NR 

REwlNil NSTIF 
REWIND NRED 
REWIND N1 
rewind N2 

MAIN l.ajP OVER ALL 8LGCKS 
OD 600 N3=lfNdEQCK * 

If- lNJ«NE«tJ GU TO 10 
READ (NSTIFI A 
IFlNLo.GTtll GD TU loO 
MAXAt 11=1 

WRITE(NRED) AtllAXA 
IFlAUll 1*1T4»3 

1 

WR1TE(6*1010 1 KiCtAlU 
3 OD u L=1 >NV 
5 A(11-LI=A(1*L|/A(11 
MR=1+NV 

WRITlCNLI tA(KK) «K=2*Ktil 
RETURN 

If (NT8*bQ.U uU TO lOtJ 
REWIND N1 
REWIND N2 
read CNl ) A 

FIND CUi-JMN HEIGHT^ 

> KU=1 

AMsMlNOIMAiNEUd) 

MAXA11)=1 
UU UO N=2*H1 
IF IN-HAI L20il2Utl30 
J KU=KU + NEDD 
K.K-AU 

KM=M1NiJ(N«KM1 
GO TO 140 


’'kZ3 r~i 'f— ! r~~i f~i ciD [— 3 ;Czii a CZD 


f 1 . 
1 

;SEsaj- 

LA5Htt. 

*FUH FORTRAN >/»oAId21) 

1 

iiOiao 

130 

KU=KU +1 


20190 


ftK=KU 

1 » 

20200 


IF iN-NcQUl i40«140»13b 


20210 

136 

HH=RH - 1 


20220 

140 

DO 160 

1 ^ 

2Q230 


IF (A(KKl) llOflbOtllU 


2024C 

160 

KK=KK - INC 

I 

20250 

ilO 

HAXA{N)=KK 


20260 

C 



20270 


IP tAdJJ 172. l74fX76 


202B0 

174 

KK-ChJ~l MfNtQJ + 1 


20290 


IF (KtC.uT.NEU) GQ Td 590 


20300 


MRITE (bilOOOt KA 

1 > 

20310 


STOP 


20320 

172 

KK=1NJ-1 l*NeOu t 1 


20330 


tfRlTi* (bilOlO) KKiAdl 


20340 

C 


[■ \ 

20 350 

C 

FACTORlZb. LEADING bi-DLA 

i '■ 

20360 

176 

OCj 200 N=2.NEQ3 


20370 


NH^HAXAIM 


20300 


IF fNri-WJ 200,200.210 

1 1 

2O390 

210 

XL=N + INC 


20400 

20410 


X=N 

0=0,00 

: 

20420 


DU 220 KK-KL*NH,1NC 

1 

20430 


K=X - i 


20440 


C=AlK)0/A(6l 


20450 


0=0 + C+AtKK) 


20460 

220 

AtKM=C 

i 

20470 


Alin=AlNl - D V 


20450 

c , 

[ 

20490 


IF (AlNd 222,224,230 

1 

20500 

224 

XX=i NJ-i J4NEQU F N 

) ■ 

20 510 


IF IKX.br.NEQJ GO TO 590 

1 

20520 


IkKl 7E 16 ,11)00 ) KK 


20530 


5TUP 


20540 

222 

KK=tNJ-l )4N£Qt) + N 

) 

2055U 


tfKlTL 10,1010} KK.,\(N} 


20560 




20570 

2j0 

lt=Ntuti 


20580 


DU 240 J=1.HA2 

1 

20590 


MJ=HAXAl NPJj - IL 

20600 


IF (HJ-N) 240*240,260 


20610 

260 

KU=MlNiUMJ,NhJ 


20620 


XNFW + it 

y . 

20630 


L=0,00 

20640 


00 300 XK=XL,Kd,^Nu 


206SO 

300 

L=C P A(XX}4AtxX-t-IC} 


20660 


AtKNI=AlKNl - C 

1 

20670 

240 

1C=IL + nLUU 


20680 - 

c 



20690 


K=N + NrtA 


2070 0 


DD 430 L=i,NV 

} 

20710 


XJ=K 

20720 


L=o,oo 


20730 


OU 44 0 aN-Xl,NH, INl. 
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t StiSOL LASHtt-.f-DR FURTRaN > /lU 2j-MAf-iTd idiots 

I '■ '■ i- ' 

20740 ' ' ^ IVJ-KJ 1 : 

20750 440 -J- A(KKJ4«<k;J) 

I 20760 At»C)-AtKI - C 

f 20770 4JO K=ft + Ntflll 

( 20780 C 

! 20790 200 CijUTINJE 

; 20800 c 

E 20810 t CARRY OVfcR INTO THaILINC ULDCKS 

2J820 OL) 400 NA=1 fHr«d 

r 20820 IF 1 1 NK+nJJ .uT.NoCUCiS.) Gu 7C-»03 

f 20840 Ut-RI 

' 20aS0 IF KNJ.8Q.1 |•OR.(NA«t£0<N^oU «t=rtSTIr 

E 20660 REAO (Nil <j 

20870 ML=HA4»E0b + I 

20880 rtR«Hiri0( lf^K.+U iNe0o*MI J 

f 20890 IF(rtA.EO.l) ML=«R 

; 20900 MDsHX - ML 

Ee 20910 KCsNcOEi <■ t N<- 1 J 4NtlJ84N6ad 

i 20920 ff=l 

209J0 C 

209*0 dU bOC M=Hi-tHR' 

20950 NH=HAAAIMJ 

r 20960 jCC-JiL + NEOa 

20970 IF IWH-Ai-J S0b,bl0»5l0 

I 20980 SIO A»NECI8 

f 20990 0=0.00 

f 21000 DQ 520 AA=Kl.fNH«lHC 

21010 L=ACI^K)/ACAI 

2102C 0=i> + C*AIAA> 

21030 AIAK>-=C 

21040 320 H=K - 1 

210SO u 1NJ=8CNJ - l> 

21060 IF IHO) B80«580*b30 

21070 530 IC=NEQtt 

,21080 OO 540 J=1,H0 

21090, MJ=MAxA(M>JJ - IC 

*: 21100 IF IMJ-AL.J b<tC,5S0.abC 

: 21110 550 KU=MlNO IHJ.NH) 

21120 ^N=S + IC 

. 21130 C=O.£J0 

21140 DU 570 KA-AE.KiJf IRL 

2S1S0 57S i,=C + At KM«A<KKFlCj 

21160 U(NN>=U{r,M] - t- 

21170 540 |C=1C F rtkUU 

21180 u 

21190 580 AN=6 t N*A 

21200 F.=NLOtl f NttA 

21210 DU oiO L=l«r<V 

21220 Aj-A 

21230 C-O.DO 

21240 DU o20 KA-AL,«r<H« iric 

21250 C=C F AIAAIAAIAJI 

21260 o20 KJ=AJ ”1 

21270 LJ{KRi=UtAN) - 

21200 NN-KA + NEUU 

21290 ‘uID K=tC F NC.UR 
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*FDK FUJ^TRAN Vj-bAt62l> /Nl *:j-MAr-7o 


Et.JOO 

c . 



ai3L0 

505 

HO.’sMO - 1 


eiaao 

500 

N=H + 1 


ai33t: 

C 


J 

21340 


IF (NT6*Nb. 11 60 TLi 560 


21350 


ftFilTb ImHED) A)HAXA 


21360 


OU 570 l=liNAV 


21370 

570 

Ali;=ti(IJ 


21 JBU 


uL TO 600 


2139U 

560 

HFtlT6(t'12J u 


21400 

C 



21410 

400 

CuMT INliE 


21420 

L 



21430 


M=Mi 


21440 


N1=N2 


21450 


N2=W 


21460 

590 

WRITE tMHtU) AtMAAA 


21470 

C ' 



214S0 

500 

CD.HTINUb 


21490 

L 



21500 

i. 

VEbTOR 3ACK5U6ST ITUTIDN 


2151 U 


□ 0 700 K-l*Nt(VV 


21520 

700 

Q(h)'=^OaDO 


21530 


REWIND NL 

t 

21540 

C 



215bO 


Ou 600 NJ=liNBt-QCK 


21560 


6Ai-K3PAi.E HKED 


21570 


REWIND NRED 


21560 


N ju^^NULlJCK-^NOF 1 


21590 


OU 601 NJA-laNJd 


21600 

dOl 

read (NHc-Dl AiHaaa 


21610 


backspace nheo 

• 

2162 0 


a=nebt 


21630 


DO 610 (_=1«NV 


21640 


DU 620 L^tNhd 


2 1 650 


6iKl=BlR-NE(lB) 

1 

21660 

820 

K=R - 1 


21670 

310 

K=K + NcoT + NEB 


21660 


KN=0 


21690 


KK=NidA 

1 

21700 


NDl F=NEUn 


airio 


IF IHJ.Ed.ll NDlF=NEUti - 1N6LDCK.4N6DO - NtDl 


21720 


DU 655 l_=l«NV 


21730 


DU 6b0 K=1.N01F 

) 

21740 

650 

al KNFK)=A(M<FK)/AtK) 


21750 


AK=KK -»■ NED6 


21760 

66b 

KN-KN + NEDT 


21770 


lF(MA*EQ>n C*U TU 91o 

1 

21760 


ML=NEQu + 1 


21790 


Kl=NEU11 


21600 


dU 6a0 H-=M|_*M1 


21610 


R.i.= Kjl- *' NtDtl 

1 

21820 


NU-MAAA(M) 


21830 


IF CRU-RLJ U60«670f670 


21640 

670 

K-NE(Jb 


21860 


AH=R 


J 




[rsz] 1==] 











ETj . =3 ■■ S ' ^ CII3 ESZj 

i\ ' " ■ 


EI3 1Z~3 ' CI^ [Til] tzzi ; t~! E:I3 CZI] 


SbSDL 

lasiiel 

•FOR FuRTRAN V.5A1021J /KL 

2d— M aT— 7u 12S.AC3 

21360 


DD doO t,=l«.N'kr 


£1870 


KJ=K 


£1800 


DO 890 KR-AUtKO* INC. 


£1890 


2{KJ)=8CxO] ^ A(XiU*cHKH) 


J £1900 

890 

KJ=RJ - 1 


£1910 


KH=KH + NEbT 


£1920 

880 

X=K ♦ NEuT 


21930 

860 

CClMTlNUe 


£1940 


N=NE08 


£1950 


!>□ 910 i-£«NEUu 


aiQbO 


Rl-=N + InC 


21970 


KU=HAXAIN} 

• 

1 £1980 


IF IKU-Kt.) 910.920*920 


2199C 

9£Q 

K=N 


22000 


UQ 930 l-=l.NV 


£2010 


KJ=K 


I £2020 


00 940 XK=N{-(k<.UilNC 


££030 


KJ=KJ - 1 


tf 22040' 

940 

8(f<lJ!=dtKJJ " A{XXJ«b(K) 


22050 

930 

K-X *- NE8T 


( £2060 

910 

N=N - 1 


££070 

c . ■■ 



22080 

91d 

KK=0 


£2090 


KN=^0. 


) 23100 


00 950 L.=liNV 


£2110 


DO 960 X=ltNE08 


£2120 


KK=XK +1 


22130 

9bO 

A(KK)=UIXN1-K) 


' 22140 

950 

NN-KN + NEdT 


22150 

C 



££160 


MRITE {.Hl. 1 (Aiiv) «K=1 iMWVl 


£2170 

800 

CONTINUE 


k 22180 

c 



22190 

1000 

FURplAT 1// 46H STOP +** ZErtu 

D1a80nAi- ENCuuNT t^hEU DDK 

2£20O 

A 

18H cOUATIOH SDlUTIUN, 

/ 

£2210 

X 

IjX*18H equation NUM6ER =* 

lo J 

F 22£20 

1010 

FORMAT 1/ 50H «ARN1Nu 4*4 NEOftTlVt: dlAOJNAI. UNCUUNIcREu 

22230 

X 

18H EQUhTIUN SOt-UriON* 

/ 

22240 

X 

13A.18H EQUATION NJHdER 

Id. dX. TNVAEUE =• D2d.< 

22250 

c 



J 22260 


RETURN 


22270 


END 



r’Auu 










} SUDPKOGI^AMb v.ALLt.3 

MlNO. 


SCA1.AKS AND AHRAVa i. NU iZAPLlClT JEFINITIUM ^ ••*’• NUT Hi-f LNefSCfc^J J. 


«S00£0 

1 

♦AN 

£ 

4H J 

J 

*NH 

4 


D 

♦ NEQ8 

o 

.SD0J7 

7 

+ MD 

lu 

• 50036 

11 

♦NRED 

1 .^ 

•S0035 

Id 

♦ Nl 

14 

4N0L.0CK 

lo 

• J 0 j4 

lu 


17 

♦ NSTIf 

£0 

.SOUJd 

£1 

>Nt(jT 

£2 

• 5 00 J 2 



2« 

*« 

£5 

•80031 

£6 

* .^OUdO 

£7 

♦K.M 

JO 


al 





jiy mniwii 









I 


5ESQL 

J-ASrii. 

..hCR 

FORTRAN 

V.5A162U /Kl 

25-MA7-76 

12:^6 

PAUc l-o 

KAXA 

J2 

♦ NttVV 

Jd 

^NEO j4 

iQ 35 

4mR 

36 

♦KR 

37 

•UAV 

40 

4NU 41 

•30042 42 

•5QC41 

43 

♦ KL 

44 

.50040 

45 

♦NrO 46 

♦hUlF 47 

4NWV 

60 

♦ IMC 

51 


52 

♦ J S3 

. 34 

'tf'lSt 

35 

*5000 r 

56 

.50006 

57 

4AA 60 

♦O 61 

■SOOUS 

63 

■ SOOQA^ 

64 

»NV 

65 

;soo0o 66 

4NWA 67 

A 

70 

.S0D02 

71 

.50001 

72 

4MK 7d 

•SQOOO 74 

4RA 

73 

%1C 

76 

9MAZ 

77 

♦MU 100 

. 101 

•50017 

102 

^KJ 

102 

«SQ016 

104 

.50015 105 

■ 5 001 4 10 6 

.SOOIu 

107 

=*KU 

110 

•50C12 

111 

.10003 1 12 

■SOJl 1 lid 

.10004 

114 

.50010 

1 15 

. 1000 j 

116 

4L 117 

• 10JO2 120 

.iuOoi 

1^1 


122 

.10000 

U3 

*m 124 

PNJA 125 

♦l- 

126 

.50027 

130 

.50026 

131 

.50025 132 

• 30024 IdJ 

4N2 

134 

.50022 

135 

. 50022 

136 

•50021 137 




TEMPORAKIES 




; 



■ AOIU& 

217 

.00000 

220 


J ■ . 


• 


SESQ!_ i NU OhTctTED 2 . 








MAlN» i.ASHCi.«FDn FUMTHAN yM tf.a-j4Alf-7d i^Sjo PAoE i 

^E2aS SUbKOUTiNc 'jDuEj 

222^#0 IHPi.lCIT rtLAl_*ai A-rt.<j-Zl L 

22300 C k 

22310 C CAUUS SES0l.»Pm«TJ.3TAfcSS ! 

22320 C CALLEP iiY i 

2233C -C ■ 

22340 C STATIC SiJl.UTlCN PMAbc. I 

22350 C ■ I 

4,2360 COMHOr* /uAR/ 

22370 CUHMQfi ZSCPaH/ NP (14 ) »N i •M4.#>V3 OT >N4.«4 

22330 CLtMHJN /oOL / Nou2v.'v»^c..;ti*LL«NFf 1F1 l|l(7 ) 

22390 C 

22400 RfcAL TTt4I ! 

22410 C 

22420 C SCiJ-yt FOR THE UISPCaOEHIENT VEt^TdRS 

22430 





ca' [~3 


C~3 a E3 




SGLEQ UAaUcU-tFGri FUflTHAN V.uAta^lJ /’Al IdS Jib PAOt *. 


H24A0 C 

22460 NS*io=.WE:iaij*Ut.^t2+(HtiANiJ-2)/N£att> 

2247(1 iFtN33a*i-T,<'isa) N3uu=naa 

224B0 N4=rti+NSd 

22490 N4=(<3FN5do 

22500 Ml = WdANJ + N2Ua -I 

t 22510 CA1.U SCSdi- tAlNLJ tAlNdl •AlN4y*N2C<9Mi3ftN0 •LMNuL.0CKfHE0b>N5u«Hi a 

[ 22520 A 4»d*2i7) i 

' 225J0 L I 

22540 C PAI-JT uIx*PUrtt-E.M5NT5 i 

22550 C 

22560 N-C=Nli-MJM^P*6 ! 

22S70 NSo=«fi&qw4HfcJl.0Ctc | 

i 22580 N4=4’42«N5a4M_ I 

[ V 22590 N5SN4-H 

f 22600 CAl-u EKHIJ(UN&) 

r 22610 CAUl. Pf([NTa(A(MU«A^t421tAlN41.1«Ncau*J40MNPiLt.tN£M.UU4i 

1 22620 A NEa«N5Df2*lt 

226-^0 C 

22640 C t-UMPoiTc ANi) PRll<lT euEMtENT SThESdcS' 

22650 C 

22660 M2-K1+4**.!. 

22670 NJi=N2+NEatJ*l_LJ 

22680 Ud=lHTOT-«31/iNbO F12J 

22690 CALI. STRCSSlAtNyi >A(N2y *A(N3)aNt£qii»L£i*i-L* i^a»NuLQLAl 

22700 L 

2271C return 

22720 END 



CUMMQN BLUC»;S 
/GAR/C+47036J 


A +0 

/ELPAB/1+501 
WP FO 

NJKNP 

♦ 1.6 

MUANU 

*17 

NCt. rvp 

+20 

N1 

+21 

N2 F<i2 

NEQ F27 

N3 

♦2j 

N4 

+24 

: 

+25 

KTUJ 

+26 

/SQL/(+l3» 
NBLOCk «-o 


♦1 

LC 

+2 

■ NP 

+4 

IF ICL 

+4 


5UBPHU0KAMS «.ALt.6D 

5XRl 55 5E5UL PrtllvTu cRRUH 

bCALArlb AND a 4RA75 i NU EXPLICIT DEFlNlllUN - »*X“ iNUl flc.FErtc.NCED X 

♦Ml I ♦NbUo 2 ♦lN3» 3 *i_o 4 All 






czd ' a c ::3 - 1=2 

^ ' : ; ■ . . ■ . ■■■ ■■■■ 


! 


SOLEO uASribL^f^C^ FUBTAAN V-:>Ali»aiJ /Kl 

t 

TEHFCRAHlES 

•A 3016 B (.qoooo & '■ .aaqui 7 

I 

SULEQ X NO ERRLR^bi JETECTEO 1 




•uojaa 20 


HAut ^-l ' 




_ . r-. r— , f 

f ,-.j I ...J L— -j f-- ■ 

1 • 


T-\ • > 

Kv^ 



1} 


) 


J 




E 






'4 f"" ' 

r j h- ' ■ vH 


STRESS .LASHEU.FljR 


FORTRAN V*SA(62iT /Ki 


2J-MAV-7U 


12: J6 PAOE t 


22730 

227AO 

22750 

22700 

22770 

22760 

22790 

22800 

22810 

22820 

22830 

228iO 

22859 

22650 

22870 

22680 

22890 

22900 

22910 

22920 

22930 

229A0 

22950 

22950 

22970 

22980 

22990 

23000 

23010 

23020 

23030 

23040 

23050 

23060 

23070 

23080 

23090 

23100 

23110 

23120 

23130 

23140 

23150 

23160 

23170 

23180 

23190 

23200 

23210 

23220 

2323U 

23240 

23250 

23260 

23270 

23280 


C 

c 

c 

c 


SU13RUUT1NE STRESS 157R«0,{3*NEabi 1.8 1 LI- tNbO»N8L.iJL)U 
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AFPEin)lX c 


Analysis of Axlsymmetrlc Laminated Cylinders 
Under Axlsyimetrlc Load 


The strain energy expression for axlsynmetrlc deformations of a 
laminated cylindrical shell which has orthotropic layers with principal 
material axes which coincide with the longitudinal and circumferential 
cylinder directions Is given by Eqs. (2.4) and (2.12) as 

r- nR [(1 + R + <L + . 


du. du.., d.+y-h. du . 

■ 7 ^^ + (1 + 

dx dx R dx 


®12 R ^ dx ^ dx ^ ^ ^22 Rh^ ^ ^ 1+ d^/R ^ 

+ ®44 

where the material properties and thlcknesser, are assured to be constant 
within each layer. 


The first variation of U can now be obtained with integration by 


parts as 
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The Euler equations for the cylinder are then 
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shell for which the pertinent equations are 
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Lat the solution for u^» ^2 ^ expressed as 
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Then the substitution of Eqs- ,(C6) Into Eqs. (C5) yields 
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from which the characteristic equation for X is given by 
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The coefflcl^ts and U 2 ere given In terms of W for each value of 


by the expressions 
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Since the solutions of Eq. (C7b) can be expressed as 
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the displacements may be expressed as 
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Consider a cylinder of length L which Is equally loaded at each 
edge by radial shearing forces. Then the radial deformation Is sym- 
metrical about the center of the cylinder length so that the coeffi- 
cients A 2 and B 2 may be set equal to zero. The boundary conditions at 
X " ± L/2 to be satisfied by the remaining coefficients are 


1 12 w 


2EiiR 


2 R ' dx 3 ' 4 R ' dx 


Oh -Q 


It should be noted, however, that only one of the first two equations 
need be satisfied since the sum of the two Is automatically satisfied 
by the solution given. By manipulation, an alternate equation to be 
satisfied Is given by 


The equations for the determination of and B^ are then 

X, a - ^) A, ,lnh X i + X,(l - ^) B. 3l«h x,^-^ 





which yield 


Stresses at top and bottom of the shell wall can now be obtained 
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The distribution of transverse shearing stresses can be obtained 


(C15a) 


Butf with a given by 
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Eq. (C15a) yields 
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The conditions of zero shearing stress at the inner and outer wall 


surfaces are satisfied since T from Eq. (C15c) is identically zero 


APPENDIX D 


differential . EQUATIONS FOR PLATE AND 
SHELL ANALYSIS 

While partial differential equations for shell layer Interface 
displacements may be obtained with the utilization of the principle of 
potential energy and Integration by parts, the same equations 
may be derived through the use of the layer equilibrium equations and 
the stres*;--stralnrdlsplacement relations for the shell. The various 
equations required are given here. Consider an element of a shell 
layer (Fig. D-1), which 1s loaded on Its top and bottom surface by 
Interlayer shear and normal stresses and on Its sides by force and 
moment stress resultants. The moments are defined with respect to the 
layer middle surface. Then the stress-resultant equilibrium equations 
may be written with the use of Eqs. (2.5.3) to (2.5.7) of Ref. 29 
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From Eqs. (Dla) and (Did) for the 1th layer the following 
expressions for the Interface shear stresses may be obtained; 
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An alternate expression for may be obtained from Eqs, (Dla) and 
(Did) for the 1+1 st layer as 
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The various equations may be readily combined for a flat plate 
to yield the following simultaneous differential equations for the 
Interface displacements 
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